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ABSTRACT 


The  Cloud/Fog  Analysis  System  (CFAS)  and  the  Cloud/Fog 
Application  Routines  (CFAR)  were  applied  to  weather  data  bases  to 
determine  their  sensitivity  to  control  parameters  and  to  type,  density, 
and  distribution  of  observing  stations.  The  data  rich  region  of  south- 
eastern United  States  was  selected  and  hourly  aviation  weather  (Service 
A) , six-hourly  synoptic  (Service  C) , and  twelve-hourly  radiosonde 
(RAOB)  observations  were  collected  for  weather  scenarios  of  interest  to 
Army  aviation.  Computer  methods  were  developed  to  process  these  data 
and  convert  them  into  a form  suitable  for  CFAS.  Objective  analyses 
and  output  displays  were  generated  using  CFAS  and  CFAR,  respectively, 
on  such  weather  variables  as  sky  cover,  lowest  cloud  base,  ceiling, 
visibility,  significant  present  weather,  and  cloud  obscuration  to  pilot's 
vision  within  discrete  flight  layers.  Results,  produced  from  large 
variations  in  the  computer  control  parameters  and  density  and  distribution 
of  stations,  were  used  to  modify  the  CFAS  and  CFAR  to  correct  for 
detected  errors  and  to  fix  the  control  parameters  so  that  users  are  now 
relieved  of  that  responsibility.  This  greatly  simplifies  the  knowledge 
and  experience  required  to  execute  CFAS.  A large  number  of  color- 

coded  displays  was  generated  to  demonstrate  feasibility,  skill,  and 
detail  that  is  possible  with  an  automated  meteorological  system  in 

providing  weather  information  tailored  to  Army  user  needs. 
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I INTRODUCTION 


1.1  GENERAL 
1 . 1 . I Overview 

The  Army  has  more  aircraft  and  pilots  than  the  Air  Force,  flics 
most  missions  in  weather  sensitive  light  displacement  aircraft  at  low 
altitudes  over  short  distances  and  times,  and  yet  provides  the  least 
amount  of  meteorological  training  for  its  pilots.  It  is  not  unusual  that 
Army  pilots  have  considerable  difficulty  and  are  error  prone  in  reading 
and  properly  interpreting  weather  codes  for  aviation,  synoptic,  radiosonde, 
and  pilot  report  observations.  The  present  standard  weather  teletype 
sequences  do  not  make  the  pilot's  job  easy.  They  have  to  memorize 
the  weather  station  call  letter  and  numbers,  rummage  through  clip  boards 
to  search  for  particular  stations  and  time  changes  of  weather,  and 
attempt  to  assimilate  the  four-dimensional  distribution  of  weather  im- 
portant to  his  flight.  The  labor-saving  and  real-time  accurate  analyses 
of  weather  by  a computer  based  Automatic  Meteorological  System  (AMS) 
offers  considerable  potential  and  value  to  the  Army  Aviation  Community. 

As  a part  of  the  AMS  program,  the  Geo-Atmospherics  Corporation 
(GAC)  developed  the  Cloud/Fog  Analysis  System  (CFAS)  to  take  any  type 
and  age  of  meteorological  observation  and  perform  an  objective  analysis 


to  specify  each  weather  variable  on  a three-dimensional  array  of  grid 
points.  This  cloud-fog  data  base  (CFDB)  generated  by  CFAS  can  be  used 
individually  or  its  parts  analyzed  to  tailor  operational  products  to  a given 
application.  The  Geo- Atmos pherics  Corporation  developed  several  Cloud- 
Fog  Application  Routines  (CFAR)  to  interrogate  the  CFAS  data  base, 
extract,  analyze,  and  display  a map  of  each  chosen  parameter,  i.e.  sky 
cover,  ceiling,  visibility,  cloud  base,  cloud  top,  significant  weather, 
and  cloud  amounts  within  nine  different  layers  extending  from  the  surface 
to  the  3,000  meter  altitude. 

This  study  was  conducted  to  determine  the  sensitivity  of  the 
Cloud  Fog  Analysis  System  (CFAS)  and  the  Cloud  Fog  Application  Routines 
(CFAR)  to  type  and  distribution  of  observations,  to  select  suitable  CFAS 
control  parameters  for  routine  use,  and  to  perform  a series  of  computer 
runs  on  weather  scenarios  to  demonstrate  potential  benefits  of  weather 
automation  to  Army  aviation. 
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1.1.2  Background 


i 


The  C FAS/C  FAR  is  a computet  software  package  of  several  mter- 
rolati  d subprograms  coded  in  the  language  of  FORTRAN  V.  The  CFAS/ 
CFAR  was  designed  to  be  om  of  the  subsystems  ot  the  U.S.  Army’s 
Automatic  Meteorological  System  (AMS).  The  function  of  the  CFAS  is  to 
create  and  maintain  a cloud-log  data  base  (CFDB)  in  near  real-time  on 
a square  grid  covering  a uscr-specitied  geographical  area  relocatable 
anywhere  in  the  world.  The  OTAR  then  utilises  the  CFDB  to  generate  on 
demand,  products  useful  m Army  aviation  operations.  CFAR  products 
currently  developed  include  two-dimensional  depictions  of  areas  ot 
critical  visibilities  on  the  ground,  ceilings,  and  severe  convective 
activity . 


• I 


which 


The  data  sources  which  the  CFAS  oojectively  analyzes  and  from 
the  CFDB  is  generated  include: 

1.  Selected  elements  trom  scheduled  teletype  network  trans- 
missions ot  surface  and  upper  an  observations  such  as 
AIRWAYS,  SYNOF,  MFTAR  and  RAPIOSON  DF  coded  messag.  s. 

2.  The  three-hour  prognosis  ot  layered  cloud  cover  produced 


by  the  Air  Force  Global  Weather  Cential’s  Three  Dimensional 


3.  Elements  of  non-scheduled  and  special  weather  reports 
corresponding  to  elements  in  either  of  the  above  sources. 

The  CFDB  consists  of  the  elements  listed  in  Table  1 . 1 specified 
at  each  grid  point  in  a horizontal  window  500  km  on  a side.  The  grid 
points  are  spaced  25  km  apart. 
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TABI.r  1.1  Elements  ot  the  Cloud  Fog  Data  Base  (CFDB) 


Element 

Units 

Total  sky  cover 

00  - 100  Percent 

Height  ot  ceiling  layer 

dekameters,  AGL* 

(minus  if  a variable  ceiling) 

Prevailing  visibility  at  surface 

meters  (minus  if  variable) 

Base  height  ot  lowest  cloud 

d e k a m et  er  s , AG  L 

Top  height  of  highest  cloud 

dekameters , AGL 

Present  weather 

WWO**  code  4077 

Percent  cloud  cover  in  the  following  layers: 

0 - 49 

meters,  AGL 

50  - 99 

it  it 

100  - 199 

it  M 

200  - 299 

• i m 

300  - 599 

it  It 

000  - 999 

it  «t 

1000  - 1499 

H H 

1500  - 1999 

•1  It 

2000  - 3000 

II  II 

*AGL  - above  ground  level 

**WMO  - World  Meteorological 
Organization 

lo 


Previous  work  on  CFAS  and  CPAK  was  directed  toward  an 


engineering  solution  to  insure  "working"  models  and  computer  programs 
for  analyzing  and  applying  cloud  tog  data.  Under  these  conditions, 
refinements  were  neither  possible,  nor  required.  Also,  in  order  to  make 
the  work  effort  tractable,  a decision  was  made  to  manually  input  para- 
meters which  control  the  CFAS  and  CFAR.  Those  input  parameters  that 
are  weather  dependent  and  that  were  varied  in  this  study  are  listed  in 
Table  1.2.  They  include  such  items  as  the  time  to  begin  an  analysis, 
time  of  the  oldest  data  to  be  used,  number  of  possible  grid  squares  to 
search  for  observations  surrounding  a grid  point,  first  search  square 
size  that  will  be  interrogated  (it  no  observation  exists,  second,  third, 
etc.,  search  square  size),  distance  between  two  or  more  observations 
to  be  combined  into  one  "best  report",  and  distance  and  time  scale 
factor  to  allow  cellular  convective  cloud  information  to  be  extravolau  1 
differently  in  space  and  time  than  the  more  strata  form  middle  and  high 
clouds. 


IT 


TABLE  1 . I 
Parameter 
TIME 
TV  MOLD 

Nb’SQ 

ISSQ(l) 

isso(-i) 

ISSQ(3) 

ISSQ(4) 

issq(s) 

DSP 

DIST(l) 

DIST(J) 

131  ST  (3) 
TYMC(l) 

TYMC(2) 

TVM  C(3) 


Weather  Related  t'l'AS  Input  Parameters 


Description 


Time  of  analysts,  minutes 

Tune  ot  oldest  data  to  be  used  in  analysis, 
minutes 

Number  ot  search  square  s 

M/e  ot  t irst  soared  ^mtan  , i^rio  units 

Si/e  ot  second  search  squute,  qua  units 

Size  ot  thud  search  square,  qua  units 

Size  ot  tourtn  search  square,  qrid  units 

Si/e  of  tilth  search  squan  , qnd  units 

Maximum  distance  between  observations  wmch  aim 
combined  into  a best  report,  Km 


Di; 

stanci 

scale 

factor 

when  e 

>nly  conveo 

t IVt 

' eh 

an 

• pm  o 

■lit  in 

a best 

report  , 

Kin 

Dr 

s taiux 

sea  li 

lactor 

W'lU  1 ! c 

Ull;  v.  Z : 1 V l C 

i i \ , 

am 

m i. 

dale  c 

lends  . 

an  pre 

‘sent  m 

>1  best  imp 

Vil  , 

Km 

Di 

stance 

scale 

tactor 

for  all 

other  eaSi 

^ , 

Km 

Time  scale  factor  wm  n o.Tt\  convective  tiouds  are 
present  in  a best  report,  minutes 

Time  soak  factor  wiu  n only  convective  and  miaale 
clouds  are  present  m a best  report,  minutes 

Time  scale  factor  tor  all  ottier  cases,  minutes 
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A 500  km  square  geographical  area  within  the  southeast  portion 
of  the  United  States  was  selected  for  analysis  and  demonstration  of  the 
C FAS  C FAR.  Three  reasons  for  selecting  this  area  are  (1)  high  data 
density,  (2)  variety  of  weather,  and  (3)  location  Army  aviation  school  at 
Fort  Rucker,  Alabama.  Attempts  were  made  to  center  the  grid  window 
symetncally  about  the  Army  aviation  training  bases  in  Alabama.  All 
routine  and  aperiodic  weather  reporting  stations  were  identified  within 
the  500  km  analysis  window  region  and  for  an  additional  200  km  border 
surrounding  that  window  to  insure  data  distribution  in  a manner  that 
allows  analyses  to  be  performed  on  the  outer  boundaries  of  the  window 
region. 

A weather  watch  was  instituted  to  search  for  and  select  synoptic 
meteorological  conditions  that  provide  sequences  and  scenarios  important 
to  Army  aviation.  Situations  depicting  deteriorating  flight  weather  con- 
ditions received  priority.  Particular  attention  was  given  to  those  features 
effecting  sky  cover,  ceiling,  visibility,  and  significant  weather  events. 

Three  different  synoptic  weather  sets  were  selected  tor  the  dates 
December  6 to  8,  1976,  February  25  to  27,  1977,  and  March  5 to  7,  1977. 
Data  included  surface  observations  of  hourly  aviation  weather  and  .six- 
hourly  synoptic  data,  and  upper  air  observations  from  radiosonde  and 
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pilot  rcjKirts . Each  data  sot  was  examined  closely  to  eliminate  errors, 
minimize  missing  data,  and  maximize  the  number  of  data  sources  and 
inputs.  Care  was  exercised  in  decoding  the  weather  data  and  encoding 
it  in  a term  compatible  with  the  exact  input  needs  of  the  CFAS/CFAK 
programs.  Inis  data  sample  was  chosen  to  provide  a technical  challenge 
to  the  C FAS/C  FAR  system  while  at  the  same  time  depicting  features  that 
are  operationally  interesting  to  pilots.  Several  types  and  changes  of 
weather  scenarios  were  captured  to  provide  a large  total  data  sample 
which  illustrates  a wide  range  of  operational  flight  weather  conditions. 

Operational  performance  and  quality  of  information  output  from 
ClAS/CTAR  was  studied  to  determine  sensitivity  to  type  and  density  of 
data  sources.  lest  runs  were  made  using  maximum  amount  of  all  types 
of  data  and  using  every  available  reporting  station.  Results  from  these 
runs  formed  the  baseline  analysis  or  ground  truth  for  subsequent  analyses 
where  data  amount  and  distribution  were  systematically  varied.  A com- 
parative analysis  was  made  to  determine  degradation  as  a function  of 
increasing  sparsity  of  data.  The  data  sparsity  was  created  in  two  ways. 
The  first  way  consisted  of  withholding  data  from  selected  stations  in 
such  a manner  as  to  uniformly  reduce  the  density  throughout  the  entire 
region.  The  second  way  consisted  of  withholding  data  from  stations  in 
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such  a manner  as  to  create  a subregion  of  little  or  no  data  (silent  area) 
and  thus  create  a nonuniform  distribution  of  data  density.  This  latter 
case  having  a silent  area  is  typical  of  battlefield  environments. 

The  present  manually  inputted  control  parameters  of  CFAS  were 
altered  to  effect  such  things  as  the  mesoscale  and  synoptic  scale  dis- 
tance and  time  constant  parameters,  cut-off  distances  for  best  reports, 
and  number  of  search  squares  allowable  to  find  data  necessary  to 
analyze  and  depict  a weather  value  at  a grid  point.  Three  methods  were 
imposed  to  establish  a quantitative  measure  of  "goodness  of  fit"  between 
'ground  truth"  and  analyses  by  sequentially  changing  control  parameters. 
The  correlation  coefficient  and  root  mean  square  error  (RMSE)  were 
computed  where  the  ensemble  included  all  values  at  every  grid  point 
where  an  analysis  was  possible.  Those  grid  points  where  an  analysis 
was  not  possible  were  tabulated  and  a percentage  of  missing  analyses 
was  computed.  A total  of  16  test  runs  was  made  using  all  available 
weather  observations  while  only  one  control  parameter  was  varied  at  a 
time.  Three  different  ranges  of  values  were  used  for  seven  distance 
and  time  control  parameters  to  bracket  the  expected  mean  with  high  and 
low  values.  The  five  search  square  control  parameters  were  varied  to 
investigate  performance  versus  computer  execution  time  required  for  an 
analysis.  A fixed  predetermined  set  of  instructions  was  defined  for 


operating  the  CFAS  and  was  incorporated  in  the  final  operating  computer 


program. 

A time  series  of  observations  was  selected  to  depict  interesting 
Army  aviation  weathe'  scenarios  for  the  Alabama  area.  The  sequence  of 
events  proceeds  from  clear  skies  to  rapidly  deteriorating  flight  weather, 
such  as  thunderstorms,  low  clouds,  and  visibility  restricting  fog  moving 
across  the  region.  The  C FAS/ C FAR  computer  analyses  show  the  rapid 
changes  that  occur  from  one  observation  time  to  another,  the  detail  of 
weather  depiction  that  is  possible,  and  the  principle  of  .displaying 

i 

weather  according  to  weather  categories  indicative  of  weather  threat  to 
flight  operations.  A packet  containing  highlights  of  these  results  was 
put  together  for  the  purpose  of  demonstrating  some  of  the  advantages 
and  value  of  automated  weather  analyses  to  Army  operations. 

; 
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2 WEATHER  DATA  BASES 


2.1  DATA  AND  ANALYSIS  WINDOWS 

A square  window  900  kin  on  a side  in  the  Southeastern  United 
States  was  selected  as  the  source  of  data  for  the  CFAS  sensitivity 
analysis,  as  shown  in  Fig.  2.1.  While  the  CFAS  program  will  accept 
data  from  any  point  within  this  900  km  window,  actual  grid  point  analyse 
occur  only  within  the  inner  500  kilometer  square  whose  southwestern 
corner  is  31°  north  latitude,  89°  west  longitude.  Figure  2.2  is  a topo- 
graphic map  of  the  500  km  inner  square  where  intersecting  lines 
represent  the  grid  points  of  points  of  analyses.  The  state  of  Alabama  is 
enclosed  within  and  encompasses  nearly  all  of  the  500  km  inner  square 
where  CFAS  performs  the  grid  point  analyses. 

These  data  windows  were  chosen  due  to  the  variability  of 
weather  conditions,  or  significant  weathers,  the  high  density  oi  weather 
service  stations  which  report  regularly  and  the  fact  that  the  Tort  Rucker 
Army  Aviation  Training  Center  is  located  in  the  center  of  the  square. 


2.2  SERVICE  A - HOURLY  WEATHER  REPORTS 

2.2.1  St.it  i « > 1 1 Location  and  Chat  actenst  ics 

I’here  are  *15  weather  service  stations  within  the  900  kilometer 
window  which  transmit  hourly  (SA)  or  supplementary  (SW)  aviation  weathei 
observations  over  teletype  service  A.  I'orty-two  stations  report  regularly 
while  3 stations  report  very  infrequently.  These  -IS  stations  report  from 
nine  states:  Alabama,  Arkansas,  Honda,  Georgia,  Louisiana,  Mississippi, 
North  Carolina , South  Carolina,  and  Tennessee.  The  station  locations 
range  from  29°  to  36°  Nortli  Latitude  and  81°  to  90°  West  Longitude. 

The  altitudes  of  the  stations  range  from  0 to  805  meters  above  sea  level. 
Table  2.1  gives  a detailed  account  of  each  Service  A station,  its 
identifier,  location,  latitude,  longitude  and  altitude  in  meters,  figure 

2.3  shows  all  Service  A stations  and  their  locations  within  the  900  kilo- 
meter window,  plotted  on  a mercator  map  projection,  1:7,  500,000  scale 
lit  latitude  22°  30".  Notice  that  the  Seivice  A stations  are  umfoimly 
distributed  throughout  the  entire  window  icgion. 


TABLE  2.1  Service  A Station  Information 


CALL 

LATITU  DE 

LONGITUDE 

ALT1TI 

LETTERS 

NAME 

NORTH 

WEST 

(METE 

ABY 

ALBANY  GA 

31° 

32' 

084° 

1 1’ 

60 

AGS 

AUGUSTA  C.A 

33° 

22' 

081° 

58’ 

45 

AHN 

ATHENS  GA 

33° 

57' 

083° 

19* 

247 

AMG 

ALMA  GA 

31° 

32' 

082° 

31' 

63 

ANB 

ANNISTON  AL 

33° 

35* 

085° 

51' 

188 

AND 

ANDERSON  SC 

34° 

30' 

082° 

43' 

236 

AQQ 

APALACHICOLA  PL 

29° 

44' 

084° 

59' 

1 1 

ATL 

ATLANTA  GA 

33° 

39' 

084° 

26’ 

315 

AVL 

ASHEVILLE  NC 

35° 

26' 

082° 

32' 

661 

BHM 

BIRMINGHAM  AL 

33° 

34' 

086° 

45' 

192 

BN  A 

NASHVILLE  TN 

3o° 

07' 

088° 

41' 

184 

BTR 

BATON  ROUGE  LA 

30° 

32' 

09  1 0 

09' 

23 

BVE 

BOOTHVILLE  LA 

2gO 

20' 

089° 

24* 

0 

CEW 

CRESTVIEW  TL 

30° 

47’ 

088° 

31' 

56 

CHA 

CHATTANOOGA  TN 

35° 

02' 

085° 

12’ 

210 

CSC. 

COLUMBUS  GA 

32° 

31' 

084° 

56' 

l 20 

CSV 

CROSSVILLE  TN 

35° 

57 1 

085° 

05' 

570 

DHN 

DOTHAN  AL 

3 1 0 

19' 

085° 

27' 

113 

DYR 

DEYERSBURG  TN 

36° 

or 

089° 

24' 

105 

FTY 

ATLANTA  C.A 

33° 

47' 

084° 

31' 

257 

GLH 

GREENVILLE  MS 

33° 

29' 

09  0° 

59' 

40 

GNV 

C.AINSVILLE  EL 

2 90 

42' 

082° 

16' 

50 

GSP 

G REEN V ILLE/S PRTN  BG  SC 

34° 

54' 

082° 

13’ 

296 

G WO 

GREENWOOD  MS 

33° 

30' 

090° 

12' 

41 

HSV 

HUNTSVILLE  AL 

34° 

39' 

088° 

46’ 

1 96 

JAN 

JACKSON  MS 

32° 

19' 

09  0° 

05’ 

101 

JBR 

JONESBORO  AR 

35° 

50' 

090° 

39' 

805 

MCB 

McCOMB  MS 

31° 

1 (>' 

09  0° 

28' 

143 

MCN 

MACON  C.A 

32° 

42' 

083° 

39' 

1 1 0 

MEI 

MERIDIAN  MS 

32° 

20’ 

088° 

45’ 

94 

MEM 

MEMPHIS  TN 

35° 

03' 

089° 

59' 

87 

MGM 

MONTGOMERY  AL 

32° 

18’ 

088° 

24’ 

62 

HBI 


TABLL’  2.1  (Continued) 


CALL 

LATITUDE 

LONG  ITU  DE 

ALTITUDE 

LETTERS 

NAME 

NORTH 

WEST 

(METERS) 

MGR 

MOULTRIE  GA 

31° 

OS' 

083° 

48' 

88 

MKL 

JACKSON  TN 

35° 

36' 

088° 

55’ 

129 

MOB 

MOBILE  AL 

30° 

41* 

088° 

15’ 

67 

MSL 

MUSCLE  SHOALS  AL 

34° 

45' 

087° 

37' 

171 

MSY 

NEW  ORLEANS  LA 

29° 

59* 

09  0° 

15’ 

9 

* NEW 

NEW  ORLEANS  LA 

30° 

02’ 

090° 

02' 

3 

PNS 

PENSACOLA  PL 

30° 

28' 

087° 

12’ 

36 

* RMG 

ROME  G A 

34° 

21 ' 

085° 

10' 

196 

* SPA 

SPARTANBURG  SC 

34° 

55' 

081° 

57' 

251 

TCL 

TUSCALOOSA  AL 

33° 

14’ 

087° 

37' 

57 

TLH 

TALLAHASSEE  FL 

30° 

23* 

084° 

22' 

21 

TYS 

KNOXVILLE  TN 

35° 

49’ 

083° 

59' 

299 

VLD 

VALDOSTA  GA 

30° 

47' 

083° 

17* 

66 

* Stations  report  never  or  only  once  in  our  data  sample 


2.2.2  Aviation  Weather  Code  and  Data  Used 


Throe  types  of  Service  A reports  wore  used  to  build  the  data 
sets.  They  were  the  regularly  scheduled  hourly  reports,  SA,  the  irregular 
reports  from  supplementary  stations,  SW,  and  the  specuil  reports,  51', 
made  when  warranted  by  significant  weather  changes.  The  format  of  each 
of  these  reports  provides  foi  the  following  information: 

1.  Station  Identifier 

2.  Type  and  Time  of  Report 

3.  Sky  Condition  and  Ceiling 

4.  Visibility 

5.  Weather  and  Obstruction  to  Vision 

6.  Sea  Level  Pressure 

7.  Temperature  and  Dew  Point 

8 . Wind  Direction  and  Speed 

9 . A1 1 1 m et er  Set t i ng 

10.  Runway  Visual  Range  or  Remarks 
For  the  purposes  of  the  sensitivity  study  we  weie  concerned  only  witii 
the  station,  type  and  time  of  report,  sky  condition  and  ceiling,  visibility 
ana  present  weather  for  those  elements  of  the  Service  A report  having 
some  bearing  on  Army  Aviation  operations,  figure  2.4  is  an  example  of 
the  actual  Service  A teletype  transmission,  while  Figure  2.5  is  a sample 

e separated  by  i ,oi 'nation  groups. 

;u> 


of  the  coded  messag 


Figure  2.4  Service  A Teletype  Transmitted  Report 
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Figure  2. 

5 Coded 

Service  A Transmission 

Location 
Identifier 
and  Type 
of  Report 

Sky  and 
Celling 

Visibility, 
Weather,  X 
Obstruction 
to  Vision 

Sea 

Level 

Pressure 

Temperature 

and 

Dew  Point 

Wi  nd 

Altimeter 

Setting 

Remarks 

CEW  SP 

-X  E3  OVC  1/2  RF 

123 

/ 65/62 

/200b 

/990/ 

RB  30 

The 

CFAS  program 

searches 

for  the  most 

current 

and  signifi 

cant 

information,  making  the  time,  information  content , and  type  of  report 
important.  A special  report  indicated  by  the  letters  SP  after  the  station 
designator  as  in  Figure  2.5,  is  weighed  more  heavily  by  CFAS  than  a 
regular  hourly  Service  A report  because  it  is  the  most  current  report  and 
because  Specials  occur  only  when  a significant  change  in  weather  has 
occurred.  In  the  Service  A transmission  the  sky  condition  and  ceiling 
inter  mat  ion  follows  the  location  and  type  ot  report,  lhe  sky  couditK  n 
is  ropoited  by  one  uf  seven  sky  covei  desigiMb  ,i  n ,||Ven  in  InMr  2..:, 


TABLE  2.2  Summary  of  Sky  Cover  Designators 


Designator 

Meaning 

Spoken 

CLR 

Clear  (Less  than  0.1  sky  cover) 

CLEAR 

SCT 

Scattered  Layer 
sky  cover) 

Aloft  (0.1  through  0.6 

SCATTERED 

BKN* 

Broken  Layer  Aloft  (0.6  through  0.9 
sky  cover) 

BROKEN 

OVC* 

Overcast  Layer 

Aloft  (More  than  0.9,  or 

1 . 0 sky  cover) 

OVERCAST 

-SCT 

Thin  Scattered] 

| At  least  1/2  of  the  sky 

THIN  SCATTERED 

-BKN 

Thin  Broken  ( 

cover  aloft  is  transparent 

THIN  BROKEN 

-OVC 

Thin  Overcast  ' 

t 

[ at  and  below  the  level  of 
' the  layer  aloft 

THIN  OVERCAST 

X* 

Surface  Based  Obstruction  (All  of  sky  is 

hidden  by 

surfaced  based  phenomena) 

SKY  OBSCURED 

-Y 

Surface  Based  1 

Partial  Obscuration  (0.  1 or 

more , but 

not  all,  of  sky  is  hidden 

SKY  PARTIALLY 

by  surface 

based  phenomena) 

OBSCURED 

* Sky  condition  represented  by  this  designator  may  constitute  a ceiling 
layer 


The  height  of  the  base  of  a layer  of  clouds  precedes  the  sky 
cover  designator.  Height  is  reported  in  hundreds  of  feet  above  ground 
level.  If  the  station  designator  is  CLR,  no  height  will  be  given  since 
no  sky  cover  is  reported.  When  more  than  one  layer  is  reported,  layers 
are  in  ascending  order  of  height.  The  height  may  be  preceded  by  the 
letter  "M"  - measured  ceilings,  "E"  - estimated  ceiling,  or  "W" 
indefinite  ceiling,  meaning  that  regardless  of  the  ipethod  of  determination. 
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vertical  visibility  is  classified  as  an  inch  unit*  coiling.  lln  ceiling 
designators  are  significant  to  CPAS  as  it  searches  tor  ami  weighs  more 
heavily  the  informal  si  presn;  . a to  m most  n liable,  t'ne  measured 
height  would  theretore  be  mon  signiticaiit  that,  eitiwi  Un  vstiu.au  : >•, 
the  indefinite  ceiling. 

The  visibility  at  the  observation  site  immediately  fob.  us  ihc 
sky  condition  and  ceiling  report.  Visibility  is  the  greatest  aistance 
objects  can  be  seen  and  identified  through  at  least  180''  ot  the  honcon . 
It  is  reported  in  statute  miles  and  fractions.  Visibility  is  important  tv 
r'l'AS  as  are  cloud  layers  or  any  other  factor  affecting  the  pilots  ability 
to  see. 

Weather  and  obstructions  to  vision  when  occurring  vu  the 
station  are  reported  immediately  following  visibility,  l'ne  term  weattni 
refers  to  those  items  listed  m Fable  2.  J rather  than  to  d.e  n.civ  „c.,  ,.,i 
meaning  ot  all  atmospluric  phenomena.  UVatiiei  is  s.gnit. card  >\  i 
regard  to  visibility  as  well  as  flying  conditions , am!  tin. refer c rs  m,- 
portant  in  thv'  CPAS  sensitivity  analysis. 


V.' 
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TABLE  2.3 


Weather  and  Obstruction  to  Vision  Symbols 


A 

BD 

BN 

BS 

D 

F 

GF 

H 

IC 

IF 

IP 

IPW 

K 

L 

R 

RW 

S 

SG 

SP 

SW 

T 

T+ 

ZL 

ZR 


Hail 

Blowing  Dust 
Blowing  Sand 
Blowing  Snow 
Dust 
Fog 

Ground  Fog 
Haze 

Ice  Crystals 

Ice  Fog 

Ice  Pellets 

Ice  Pellet  Showers 

Smoke 

Drizzle 

Rain 

Rain  Showers 
Snow 

Snow  Grains 
Snow  Pellets 
Snow  Showers 
Thunderstorm 
Severe  Thunderstorm 
Freezing  Drizzle 
Freezing  Rain 


Precipitation  intensities  are  indicated  thus: 

— Very  Light;  -Light;  (no  sign)  Moderate;  + Heavy 
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2.3  SERVICE  C - SYNOPTIC  WEATHER 
2.3.1  Station  Location  and  Characteristics 

There  are  27  weather  service  stations  within  the  900  kilometer 
window  which  transmit  Service  C observations,  17  stations  report 
regularly  while  10  stations  report  very  infrequently.  All  those  stations 
reporting  Service  C also  report  Service  A and  report  from  the  same  states 
as  Service  A,  except  Arkansas.  Again  the  station  locations  range  from 
29°  to  36°  North  Latitude  and  81°  to  90°  West  Longitude.  The  altitudes 
of  the  stations  range  from  0 to  661  meters  above  sea  level.  Table  2.4 
gives  a detailed  account  of  each  Service  C station,  its  identifier, 
location,  latitude,  longitude  and  altitude  in  meters.  Figure  2.6  shows 
all  Service  C stations  and  their  locations  within  the  900  kilometer 
window.  Notice  that  in  addition  to  being  fewer  in  number,  Service  C 
stations  are  not  as  uniformly  distributed  as  the  Service  A stations. 
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TABLE  2.4  Service  C Station  Information 


STATION  CALL 

tf  LETTERS  NAME 


LATITUDE  LONGITUDE  ALTITUDE 
NORTH  WEST  (METERS) 


* 72213 

AYS 

WAYCROSS  GA 

31° 

15' 

082° 

24' 

46 

72214 

TLH 

TALLAHASSEE  EL 

30° 

23' 

084° 

22* 

21 

* 7221b 

ABY 

ALBANY  GA 

31° 

32’ 

084° 

1 1 1 

60 

72217 

MON 

MACON  GA 

32° 

42' 

083° 

39' 

1 10 

72218 

AGS 

AUGUSTA  GA 

33° 

22' 

081° 

58' 

45 

72219 

ATL 

ATLANTA  GA 

33° 

39* 

084° 

26' 

315 

* 72220 

AQQ 

APALACHICOLA  EL 

29° 

44' 

084° 

59' 

1 1 

* 72221 

vrs 

ELGIN  APB  FL 

30° 

29' 

086° 

31' 

29 

* 72222 

PNS 

PENSACOLA  TL 

30° 

28' 

087° 

12' 

3 b 

72223 

MOB 

MOBILE  AL 

30° 

41  * 

088° 

15* 

b7 

* 72224 

7222b 

MG  M 

MONTGOMERY  AL 

32° 

18' 

08  b° 

24' 

u2 

* 72227 

72228 

BH  M 

BIRMINGHAM  AL 

33° 

34' 

086° 

45' 

192 

72231 

MSY 

NEW  ORLEANS  LA 

29° 

59 1 

090° 

15' 

9 

* 72232 

BVE 

BOOTH  VILLE  LA 

29° 

20' 

089° 

24' 

0 

72234 

ME1 

MERIDIAN  MS 

32° 

20' 

088° 

45’ 

94 

72235 

JAN 

JACKSON  MS 

32° 

19* 

090° 

05' 

101 

7231  1 

AHN 

ATHENS  GA 

33° 

57' 

083° 

19* 

2 1/ 

72312 

GSP 

GREENVILLE /SPRTNBRG  SC  34° 

54' 

08  2° 

13’ 

) ■ 1 . 

* 72313 

SPA 

SPARTANBURG  SC 

34° 

55' 

081° 

57' 

2M 

723  15 

AVL 

ASHEVILLE  NC 

35° 

26' 

082° 

32’ 

66  1 

* 72323 

HSV 

HUNTSVILLE  AL 

34° 

39' 

086° 

46' 

1 9 6 

72324 

CHA 

CHATTANOOGA  TN 

35° 

02' 

085° 

12’ 

210 

72326 

TYS 

KNOXVILLE  TN 

35° 

49* 

083° 

59' 

299 

72327 

BNA 

NASHVILLE  TN 

36° 

07’ 

086° 

41' 

184 

72334 

MEM 

MEMPHIS  TN 

35° 

03’ 

089° 

59’ 

87 

* Those 

stations  not 

reporting  in  our  data 

sample 

87 


2.3.2 


S VIM  OP  Weather  Code  ami  Data  Used 


The  complete  form  of  surface  synoptic  report  used  in  the  CFAS 
sensitivity  study  is  referred  to  as  the  primary  synoptic,  the  o-hourly 
report  or  SYNOP.  The  primary  synoptic  is  ieported  at  the  standard  hours 
of  observation  which  are  0000,  0600,  1200,  1800  GMT.  The  repoits  au 
divided  into  the  universal  group  and  the  supplementary  group.  The  ttrst 
seven  or  eight  groups  of  the  SYNOP  code,  llin  through  6P0  i* , 1',  , aie 

the  universal  or  mandatory  groups.  That  information  selected  tiom  SYNOP 
for  the  sensitivity  analysis  is  from  the  universal  groups.  The  complete 
symbolic  form  ot  message  used  by  United  States  stations  in  the  -19 
continental  states  for  the  primary  synoptic  report  is  represented  m 
Fig.  2.7. 

Figure  2.7  Symbolic  Primary  Synoptic  Report 

Iliii  Nddff  WwwW  PFPTT 

NhCLhCMCH  TdTdaPP 
6Popopopo  (?RRR+s)  (8NsChshs) 

(9SpSpSpSp)  (2R24R24R24R24) 

(3PwrWHwHw)  (dwUwFwHwHw) 

(4TxTxTnTn)  (Additional  Plain  Language  Data) 


Figure  2.8  is  an  example  of  an  actual  Service  C teletype  transmitted 
report.  Only  a meteorologist  who  works  daily  with  these  codes  can 
remember  and  rapidly  interpret  significant  weather  information  contained 
within  this  string  of  numbers. 


Figure  2.8  Service  C Teletype  Transmission 
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Table  2.5  defines  those  symbols  of  the  SYNOP  code  used  in  the  CPAS 
computer  sensitivity  analysis. 


i 
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TABLE  2.5  Definition  ot  SYNOP  Code  Symbols 


Iliii 

N 


CL 

h 

Cm 

Ch 

w 

w 

wvv 


Numeric  code  for  reporting  station 
Fraction  of  celestial  dome  covered  by  clouds 
Fraction  of  celestial  dome  covered  by  all  the  low 
cloud  (s)  present.  If  no  low  cloud(s)  is 
present,  that  fraction  covered  by  all  the 
middle  cloud (s)  present 
Type  of  low  cloud  if  low  cloud  is  present 
(Sc,  St,  Cu,  C^) 

Height,  above  ground,  of  the  base  of  lowest  cloud 
seen 

Type  of  middle  cloud  if  middle  cloud  is  present 
(Ac,  As,  Ns) 

Type  of  high  cloud  if  high  cloud  is  present 

(Ci,Cc,  Cs) 

Horizontal  visibility  at  surface 
Past  weather 
Present  weather 
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Service  C reports  are  similar  to  Service  A reports  in  that  they 
arc  surface  based  observations  and  report  much  of  the  sumo  nuorin.it n-n, 
such  as,  visibility,  present  weather,  sea  level  pit  s sure,  temp,  a i a.i  o 
and  dew  point  temperature.  For  the  purposes  ot  the  sensitivity  ar.alvM.-i 
the  important  Service  C information  which  supplenu  nts  the  Service  A 
reports  is  the  detailed  cloud  information. 
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2.4  RADIOSONDE  - UPPER  AIR  DATA 


2.4.1  Station  Location  and  Characteristics 

There  are  8 weather  service  stations  within  the  9 00  kilometer 
window  which  transmit  radiosonde  observations  (KAOCs),  / stations 
reported  regularly  while  one  never  reported  in  our  data  sample.  These  7 
stations  report  from  five  states:  Florida,  Georgia,  Louisiana,  Mississipp 
and  Tennessee.  The  station  locations  range  from  29°  to  36u  North 
Latitude  and  82°  to  90°  West  Longitude.  The  altitudes  of  the  stations 
range  from  0 to  32t>  meters  above  sea  level.  I able  2.ti  gives  a detailed 
account  of  each  station  which  transmits  radiosonde  reports,  the  station 
identifier,  location,  latitude,  longitude  and  altitude  in  meters.  Figure  2. 
shows  the  stations  and  their  locations  within  the  900  kilometer  window. 
Notice  that  there  are  no  radiosonde  stations  within  the  entire  state  of 
Alabama,  which  corres ponds  rather  closely  with  the  inner  SOU  km  window 
where  all  detailed  computer  analyses  are  made. 


TABLE  2.6  Radiosonde  Station  Information 


STATION  CALL 
# LETTERS 


NAME 


LATITUDE  LONGITUDE  ALTITUDE 
NORTH  WEST  (METERS) 


72213 

AYS 

WAYCROSS  G A 

31° 

15' 

082° 

24' 

46 

72220 

AQQ 

APALACHICOLA  FL 

29° 

44 ' 

084° 

59' 

11 

72221 

VPS 

ELGIN  ArB  FL 

o 

o 

CO 

29' 

086° 

31' 

29 

* 72227 

MARIETTA/ 

DOBBINS  APB  GA 

33° 

55' 

084° 

31' 

326 

72232 

BVE 

BOOTHVILLE  LA 

29° 

20' 

089° 

24' 

0 

72235 

JAN 

JACKSON  MS 

32° 

19' 

09  0° 

05' 

101 

7231  1 

AHN 

ATHENS  GA 

33° 

57' 

083° 

19' 

247 

72327 

BNA 

NASHVILLE  TN 

36° 

07' 

086° 

41' 

184 

* Never  reports  in  our  data  sample 

2.4.2  RAOB  Weather  Code  and  Data  Used 

Upper  air  soundings  or  radiosonde  observations  (RAOB)  are  reported 
in  the  radiosonde  message.  The  standard  times  of  RAOB  messages  are 
0000,  0600,  1200,  and  1800  GMT.  The  reports  are  divided  into  four 
parts,  PARTS  A,  B,  C,  and  D,  each  one  able  to  be  handled  as  a 
separate  message.  Parts  A and  B report  only  data  for  levels  up  to  and 
including  the  100  mb  level.  Data  for  levels  above  the  100  mb  level  are 
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reported  in  Parts  C and  D.  Parts  A and  B report  the  levels  of  observa 
tion  of  Army  aviation  interest,  consequently,  for  the  purposes  of  the 
CFAS  sensitivity  study  we  were  concerned  only  with  Parts  A and  B. 

The  complete  symbolic  form  of  radiosonde  message  used  by  United  Star 


stations  for  Parts  A and  B is  represented  in  Figure 


2.10, 


figure  2.10 

Symbolic  F'orm 

of  Radiosonde  Message 

PART  A 

TTAA 

YYGGId 

Iliii 

99P0P0P0 

T0T0TaoR)  Do 

dodotofofo 

OOhhh 

TTTaDD 

ddfff 

BShhh 

TTTaDD 

ddfff 

7 Ohhh 

TTTaDD 

ddfff 

SOhhh 

TTT,DD 

ddfff 

4 Ohhh 

TTTaDD 

ddfff 

JOhhh 

TTTaDD 

ddfff 

25hhh 

TTT,DD 

ddfff 

20hhh 

TTTaDD 

ddfff 

15hhh 

TTT.DD 

ddfff 

1 Ohhh 

TTTaDD 

ddfff 

88PtPtPt 

V^atDtDt 

1-1 

dtdt£tftft 

L 

77^ 

or  1 Pm  Pm  Pm 
66  ) 

dmdm^rn^m^m 

4vbvbvava  j 

PART  B 

TTBB 

YYGG/ 

Iliii 

''OPqPqPq 

Po^o^aoDo  Dq 

11PPP 

TTTaDD 

22PTP 

TTTaDD 

33PPP 

TTTaDD 

44PPP 

TTTaDD 

51515 

l°1AdfAdf 

1 

u> 


Standard  isobaric  surfaces  are  referred  to  as  mandatory  levels 


which  are  reported  in  PART  A of  the  radiosonde  message.  Data  for  these 
mandatory  levels  used  by  CFAS  are  the  standard  isobaric  surface  indica- 
tor, geopotential  height,  temperature  and  dew  point  depression. 

Significant  levels  are  defined  as  levels  at  which  temperature 
and/or  relative  humidity  data  are  sufficiently  important  or  unusual  to 
warrant  the  attention  of  the  forecaster,  and/or  are  required  for  precise 
plotting  of  the  radiosonde  observation.  These  significant  levels  are 
reported  in  PART  B of  the  radiosonde  message.  A sufficient  number  of 
significant  levels  must  be  reported  so  that  a linear  interpolation  between 
any  two  consecutively  transmitted  levels  will  give  a close  approximation 
to  the  observed  data.  Once  again  for  the  purposes  of  the  sensitivity 
study  the  level  indicator,  temperature  and  dew  point  depression  are  the 
important  pieces  of  information.  With  this  information  a determination  as 
to  whether  or  not  a cloud  layer  exists  at  each  mandatory  or  significant 
level  can  be  made.  CFAS  will  consider  this  RAOB  information  along  with 
the  Service  A and  C information  and  weigh  that  information  presumed  to 
be  most  reliable  and  combine  that  information  unique  to  each  .creating 
the  most  accurate  fully  detailed  synoptic  picture  possible,  before  making 
the  final  grid  point  analyses.  Figure  2.11  is  an  example  of  a RAOB 
teletype  transmission.  Table  2.7  defines  those  symbols  of  the  radiosonde 


TABLE  2.7  Definition  of  Radiosonde  Code  Symbols 


lliii  Numeric  code  for  reporting  station 

99  Indicator  figures . Identifies  the  surface  data  jioups 

(i.e.  , 99P0P0P0) 


hhh 


Geopotential  in  geopotential  meters,  or  tens  of 
geopotential  meters,  of  the  standard  isobaric  surface 
specified  by  the  surface  indicator.  * 


ToTo 


Observed  temperature  of  the  air  in  whole  degrees 
Celsius  at  the  surface 


ao 


Approximate  tenths  value  and  pius  or  minus  sign 
indicator  of  the  air  temperature  at  the  surface 


00,  85,  70, 
50,  40,  30, 
25,  20,  15 
and  1 0 


Standard  isobaric  surface  indicators.  Identify  the 
data  groups  for  the  1000,  850,  700,  o00,  400,  300. 
250,  200,  150,  and  100  mb  levels  in  PART  A 


TT 


Observed  temperature  of  the  air  in  whole  degrees 
Celsius  at  the  specified  pressure  level 


Approximate  tenths  value  and  plus  or  minus  sign 
indicator  of  the  air  temperature  at  the  sp. 
isobaric  surface  or  significant  level 


1 1,  22,  33,  etc.  Significant  level  indicator  numbers  in  PART  b 


*Geopotential  heights  are  reported  in  whole  geopotential  meters  up  to 
but  not  including  the  500  mb  surface.  Geopotential  heights  are  reported 
in  tens  of  geopotential  meters  (i.e.  , decameters)  for  the  500  ml.  surface 
and  higher 


code  used  in  the  CFAS  sensitivity  analysis. 


Figure  2.11  RAOB  Teletype  Transmitted  Report  (PARTS  A and  B) 
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2.5  VARIABLE  STATION  DENSITY  DISTRIBUTION 

Several  station  density  distributions  were  used  in  the  perfor- 
mance of  the  sensitivity  analysis.  The  maximum  station  density  available 
in  the  "Alabama  Square"  consisted  of  all  of  the  Service  A Airways  station^ 
included  within  the  square.  The  fractional  density  distributions  were 
created  by  eliminating  selected  stations  from  the  data  set. 

The  several  station  densities  selected  are  listed  in  Table  2.8 
together  with  the  number  of  stations  in  the  distribution  and  the  mean  dis- 
tance between  stations.  Figures  2.12  through  2.18  illustrate  the  various 
station  distributions  listed  in  Table  2.8.  Tables  2.9  through  2.14  show 
for  each  of  the  station  distributions  a list  of  the  stations  used  and  for 
each  station  the  closest  neighboring  station  and  the  distance  between 
them. 

The  stations  selected  for  elimination  to  produce  the  fractional 
density  distributions  were  chosen  subjectively  such  that  the  resultant 
station  distribution  was  as  uniformas  possible.  Starting  with  Figure  2.12, 
the  full  station  density,  the  analyst  eliminated  1/4,  1/2,  3/4,  7/8  and 
15/16  of  the  full  density  stations  to  produce  respectively  the  3/4,  1/2, 
1/4,  1/8,  and  1/16  density  station  distributions.  In  the  case  of  the 
3/4  density  station  distribution  we  compared  the  results  of  the  subjective 
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method  with  an  objective  technique.  In  the  objective  technique  we  elimi- 
nated 13  stations  from  the  full  station  distributions  which  were  closest 
to  one  another  as  per  Table  2.9.  The  results  obtained  by  the  objective 
technique  are  shown  in  Figure  2.19  and  Table  2.15.  It  will  be  noted  that 
there  is  fairly  good  agreement  between  these  results  and  those  obtained 


by  the  subjective  method. 
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TABLE  2.8  Station  Densities  Used  in  CFAS  Sensitivity  Analysis 


Density 

No.  of  Stations 

Mean  Distance 
Between  Stations 

Full 

41 

93  km 

3/4 

27 

1 19  " 

1/2 

19 

131  " 

i/4 

10 

198  “ 

1/8 

5 

27  7 " 

1/16 

2 

540  " 

1 

S.A. 

20 

N.A. 

\ 

S.A.  - Silent  Area 
N.A.  - Not  Applicable 


TABLE  2.9 


Full  Station  Density 


Station 

Nearest 

Neiqhbor 

Distance 

Km 

Station 

Nearest 

Neiqhbor 

Distance 

Km 

ABY 

MGR 

55 

GLH 

G WO 

83 

AGS 

AHN,  AND 

150 

GNV 

VLD 

150 

AHN 

AND 

90 

GSP 

AND 

65 

AMG 

VLD 

110 

G WO 

GLH 

83 

ANB 

BHM 

85 

HSV 

MSL 

81 

AND 

GSP 

70 

JAN 

MEI 

125 

AQQ 

TLH 

96 

JBR 

MEM 

no 

ATL 

FTY 

15 

MCB 

BTR 

98 

AVL 

GSP 

68 

MCN 

CSG 

125 

BHM 

ANB 

82 

MEI 

JAN 

125 

BN  A 

CSV 

150 

MEM 

JBR 

no 

BTR 

MCB 

98 

MGM 

CSG,  DHN 

145 

BVE 

MSY 

105 

MGR 

ABY 

55 

CEW 

PNS 

71 

MKL 

DYR 

68 

CHA 

CSV 

110 

MOB 

PNS 

100 

CSG 

MCN 

125 

MSL 

HSV 

81 

CSV 

TYS 

103 

MSY 

BTR 

105 

DHN 

CEW 

115 

PNS 

CEW 

71 

DYR 

MKL 

68 

TCL 

BHM 

92 

FTY 

ATL 

15 

TLH 

AQQ 

96 

VLD 

MGR 

68 

60 


TABLE  2.10  3/4  Station  Density 


Station 

Nearest 

Neiqhbor 

Distance 

Km 

Station 

Nearest 

Neiqhbor 

Distance 

Km 

ABY 

VLD 

118 

DHN 

ABY 

121 

AHN 

ATL 

110 

HSV 

MSL 

81 

ANB 

FTY 

128 

MEI 

TCL 

150 

AMG 

VLD 

110 

MEM 

DYR 

125 

AQQ 

TLH 

96 

MCN 

CSG 

125 

ATL 

FTY 

15 

MGM 

CSG,  DHN 

140 

BNA 

HSV 

170 

MOB 

MEI 

188 

BTR 

MSY 

105 

MSY 

BTR 

105 

CHA 

TYS 

145 

TCL 

MEI 

148 

CSG 

MCN 

125 

TLH 

AOO 

96 

DYR 

MEM 

125 

TYS 

CHA 

145 

FTY 

ATL 

15 

VLD 

TLH 

no 

GSP 

AHN 

150 

MSL 

HSV 

81 

GWO 

MEM 

179 

TABLE  2 

Station 

.11  1/2  Station 

Nearest 

Neiqhbor 

Density 

Distance 

Km 

Station 

Nearest 

Neiqhbor 

Distance 

Km 

ANB 

ATL 

135 

GSP 

AVL 

70 

AMG 

TLH 

212 

GWO 

MEI 

181 

AQQ 

TLH 

96 

HSV 

MSL 

81 

ATL 

CSG, ANB 

140 

MEI 

TCL 

150 

AVL 

GSP 

70 

MGM 

CSG,  DHN 

140 

BNA 

HSV 

170 

MSL 

HSV 

81 

BTR 

MSY 

105 

MSY 

BTR 

105 

CSG 

MGM,  DHN, 

ATC 

140 

TCL 

MEI 

150 

DHN 

MGM,  CSG 

140 

TLH 

AQQ 

96 

DYR 

MSL 

222 

61 


■ 


TABLE  2.12  i/4  Station  Density 


Station 

Nearest 

Neighbor 

Distance 

Km 

Station 

Nearest 

Neiqhbor 

Distance 

Km 

AQQ 

DHN 

180 

DHN 

MGM 

140 

ATL 

MGM 

240 

GWO 

MEI 

181 

AVL 

ATL 

271 

HSV 

BNA 

170 

BNA 

HSV 

170 

MGM 

DHN 

140 

BTR 

MEI 

305 

MEI 

GWO 

181 

TABLE  2, 

.13  1/8  Station  Density 

Nearest 

Distance 

Nearest 

Distance 

Station 

Neiqhbor 

Km 

Station 

Neiqhbor 

Km 

AQQ 

MGM 

315 

HSV 

ATL 

248 

ATL 

MGM 

240 

MGM 

ATL 

240 

GWO 

HSV 

341 

TABLE  2. 

14  l/l  6 Station  Density 

Nearest 

Distance 

Nearest 

Distance 

Station 

Neiqhbor 

Km 

Station 

Neiqhbor 

Km 

ATL  G WO  540 


GWO  ATL 


540 


TABLE  2.L5  3/4  Station  Density  by  Objective  Technique 


Nearest  Distance 

Station Neighbor Km Station 


Nearest 

Neighbor 


Distance 

Km 


ABY 

DHN.VLD 

120 

CWO 

JAN 

130 

AGS 

AND 

150 

HSV 

BHM 

125 

A MG 

VLD 

110 

JAN 

MEI 

125 

AND 

AVL 

105 

JBR 

MEM 

110 

ATL 

MCN 

130 

MCN 

CSC 

125 

AVL 

AND 

105 

MCI 

JAN 

125 

BUM 

HSV 

125 

MEM 

JBR 

110 

BNA 

CSV 

150 

MGM 

CSC , DHN 

140 

BTR 

MSY 

105 

MKL 

MEM 

122 

BVE 

MSY 

105 

MOB 

PNS 

1 00 

CHA 

CSV 

110 

MSY 

BTR 

105 

CSC. 

MSN 

125 

PNS 

MOB 

100 

CVS 

CHA 

110 

TLH 

VLD 

105 

DUN 

ABY 

120 

VLD 

AMC. , TLH 

1 05 

GNV 

VLD 

150 

Thus,  Figures  2.12  to  2.17  show  graphically  the  number  and 
distribution  of  weather  stations  used  to  study  the  effect  of  data  density 
on  the  objective  analysis  results.  These  cases  are  more  representative 
of  the  type  oi  data  distributions  expected,  ranging  from  highly  to  sparsely 
populated  regions.  Figure  2.18,  on  the  other  hand,  is  more  representative 
of  a battlefield  situation  where  observations  exist  in  the  friendly  territory 
but  none  exist  from  the  opponent's  region,  called  a silent  area.  The 
influence  of  these  data  distributions  on  the  CFAS  analysis  program  is 
discussed  in  section  4. 


2.6  SELECTED  WEATHER  SEQUENCES 


This  section  contains  plots  of  weather  observed  within  the 
900  km  square  centered  about  the  state  of  Alabama.  Eight  time  sequences 
of  weather  are  plotted  beginning  26  February  1977  at  1640Z  and  ending 
27  February  1977  at  0810Z,  as  shown  in  Figures  2.20  to  2.27.  This  time 
series  was  used  in  analyzing  and  preparing  material  which  demonstrates 
capabilities  of  automated  weather  depiction  for  pilots.  Anotner  weather 
plot  is  given  for  7 March  1977  at  0026Z  since  it  represents  data  used  in 
the  sensitivity  study  on  the  effects  of  changes  in  control  parameters  and 
type,  distribution,  and  density  of  observations.  Individual  weather  plots 
were  made  at  each  station  location  because  the  standard  facsimile  surface 
weather  charts  were  not  sufficiently  legible,  which  is  often  the  case  when 
attempting  to  read  printed  detail  from  facsimile  charts.  The  weather  data 
were  plotted  from  the  teletype  paper  strips  as  follows: 

1.  scattered  and  broken  clouds  are  indicated  by  a single  ana 
double  vertical  line  through  the  station  circle, 

2.  overcast  skies  are  represented  by  a solid  circle, 

3.  visibility  is  given  in  miles  to  the  left  of  the  station  circle, 

4.  ceiling  is  given  in  hundreds  of  feet  below  the  station 
circle,  and 

present  weather  symbols  are  located  on  the  lower  left  side 

66 
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of  the  station  circle.  Three  horizontal  lines  indicate  fog,  the  figure 
eight  laying  on  its  side  implies  haze,  a comma  and  period  represent 
drizzle  and  rain,  respectively,  a triangle  depicts  showers,  and  the 
capital  letter  R with  an  arrow  represents  a thunderstorm. 

Progression  of  weather  throughout  the  February  time  sequences 
varied  considerably  from  good  to  bad  flight  conditions.  The  time  series 
begins  26  February  1977  with  1600Z  routine  hourly  weather  observations 
and  special  observations  up  to  40  minutes  past  the  hour  plotted  on  the 
same  chart.  Whenever  a special  occurred,  it  was  plotted  in  lieu  of  the 
hourly  value  to  show  the  most  current  weather.  One  of  the  reasonings 
for  presenting  data  in  this  manner  is  because  the  computer  can  rapidly 
update  an  analysis  to  account  for  significant  weather  changes.  This 
feature  of  rapid  updating  and  displaying  most  current  weather  is  a distinct 
advantage  made  possible  by  the  computer  analysis  system.  Rather  than 
have  a pilot  hover  over  the  teletype  output  searching  for  specials  ot 
interest  to  him,  the  computer  system  can  accomplish  this  task  auto- 
matically. plus  properly  assimilate,  incorporate,  and  tailor  this  informa- 
tion to  his  needs.  Thus,  many  of  the  weather  plots  were  purposely 
selected  at  times  other  than  those  corresponding  to  the  standard  on-the- 
hour  report.  Not  only  does  this  allow  for  and  incorporate  more  non- 
routine special  observations  but  it  is  actually  more  realistic  in  terms  of 

67 


what  pilots  experience.  That  is  to  say,  a pilot  may  prepare  for  his 
flight  at  any  time  during  the  day  and  it  would  be  unusual  to  correspond 
precisely  with  the  fixed  time  of  a routine  meteorological  observation. 

Notice  that  the  2(>  February  1640Z  map  begins  the  weather 
sequence  with  clear  skies  in  the  eastern  and  cloudy  conditions  in  the 
western  section  of  map.  Visibilities  are  high  at  all  locations  with  the 
only  reixirted  restriction  to  visibility  being  haze  at  one  location.  Ceilings 
range  from  2000  feet  to  unlimited  in  Alabama.  Tour  hours  later  the  first 
rain  showers  and  thunderstorms  begin  to  appear  at  the  westernmost 
stations,  yet  the  visibilities  remain  high.  Three  hours  later  at  0000Z 
February  27,  showers  and  thunderstorms  begin  to  restrict  visibility,  lower 
ceilings,  and  cause  fog  resulting  from  precipitation.  A pattern  becomes 
more  discernible  at  0340Z  February  27  showing  a northeast-southwest 
axis  of  rain,  showers,  and  thunderstorms.  This  pattern  prevails  through  - 
out the  remaining  time  sequence,  moves  slowly  eastward,  and  inlensifbs. 

Ahead  of  the  squall  line  m the  clear  eastern  section,  radiation 
fog  is  forming  and  restricting  visibilities  under  clear  skies.  The  radiation 
fog  is  pronounced  only  in  the  southeastern  portion  of  the  map  and  reduces 
visibilities  to  as  low  as  one  quarter  of  a mile.  Many  changes  in  meso- 
scale  weather  features  are  contained  within  this  time  sequence.  Details 


and  changes  will  become  more  apparent  when  viewing  results  in  chapter 
5 provided  by  computer  analyses. 

The  map  for  7 March  1977  0026Z,  Figure  2.28,  constitutes  the 
verification  data  used  to  study  sensitivity  of  analyses  due  to  type, 
distribution,  and  density  of  weather  data.  In  addition  to  special  weather 
observations  at  0026Z,  the  available  computer  data  base  included  all 
standard  and  non-routine  observations  for  the  previous  12  hours.  The 
horizontal  distribution  of  weather  shows  clear  skies  in  the  northwest 
section  of  the  map,  cloudy  conditions  along  a northeast-southwest  line, 
and  haze,  fog,  drizzle,  and  rain  dominating  the  entire  southeastern 
sector.  Visibilities  ranged  from  one  mile  in  rain  and  fog  to  15  miles 
under  clear  skies.  Ceilings  ranged  from  being  non-existent  during  clear 
and  partly  cloudy  skies  to  a low  of  300  feet  in  fog  and  rain.  Thus, 
weather  events  depicted  on  the  verification  map  show  considerable 
diversity  in  type  and  range  of  values. 
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2.7  DATA  ENTRY  PROCEDURE 


The  Service  A aviation  hourlies,  Service  C Synoptic  and  Upper 
Air  data  were  obtained  in  the  medium  of  teletype  printed  paper  from  the 
National  Weather  Service  Field  Office  in  Atlanta,  GA.  A problem  that 
had  to  be  solved  was  that  of  putting  a rather  large  amount  of  data  onto 
a computer  card  or  magnetic  tape  in  the  format  required  by  the  CPAS. 

Our  approach  was  to  enter  the  data  into  the  buffer  memory  of 
our  Decwriter  II  LA-3t>  terminal  in  a format  convenient  to  the  operator. 
'Hus  enabled  the  data  to  be  entered  in  considerably  less  time  than 
would  have  been  required  had  we  used  the  CPAS  format.  The  Decwi iter's 
memory  was  sufficient  to  allow  several  hours  worth  of  Service  A,  r and 
Upper  Air  data  to  be  entered.  However,  for  convenience,  we  entered  onl . 
an  hour's  worth  at  a time.  The  off-line  text  editing  capability  of  the 
Decwriter  was  used  to  check  and  correct  eirors  in  the  data  entry.  Alter 
checking  and  correcting  any  errors,  each  hour's  worth  of  data  vva  , out- 
put to  a cassette  tape  and  the  memory  cleared  for  entering  another 
hour's  worth  ot  data.  The  format  used  in  entering  the  three  types  of 
data  is  described  in  Tables  2.16,  2.17,  and  2.18. 

Several  hours'  worth  of  data  were  accumulated  on  a cassette. 

The  entire  data  sample  required  fourteen  cassettes.  The  cassette  data 


were  input  to  mass  storage  files  on  a Decsystem  10  computer.  A 
FORTRAN  10  program  was  written  and  implemented  on  the  Decsystem  10 
to  convert  the  data  from  the  fast  entry  format  to  the  format  required  by 
the  CFAS.  The  CFAS  formatted  data  were  then  output  to  nine  track 
binary  tapes  which  could  be  read  by  the  Univac  110b  on  which  our 
version  of  the  CFAS  is  implemented. 

Using  the  procedure  outlined  above  a data  base  consisting  of 
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TABLE  2.16  Format  for 

1 

Fast  Entry  of  Service  A Airways  Da  .a 

= 

Line 

Data  List 

4 

1 

MMdd,  HHmni 

jj 

jj 

2 

III.  Ty,  Nc , Nw,  V 

3 to  3 + N^-l 

hhh,  N s , hc 

• 

, 

3 + Nc 

wlf  w2 w8 

3 +•  Nc  +•  l 

PPP,  TT,  TdTd 

* 

i 

MMdd  - Four  digit  date  code  MM  - Month  of  year,  01-12 

dd  = Day  of  month,  01-31 

HHmni  - four  digit  time  code  HH  = Hour  of  day,  00-24 

mm  - Minute  of  hour,  00  Ou 

III  - Three  letter  station  identifier 

1\  Type  code  A = SA,  hourly  airways 
S = ST,  special 

Nc  - Number  of  cloud  layers  reported,  enter  1 if  a clear  sky 
Nw  - Number  of  weathers 

V - Visibility  in  miles,  fractions  coded  as  a floating  point  decimal 

hhh  - Height  of  base  of  cloud  layer  in  hundreds  of  feet 

Ns  - Literal  description  of  cloud  cover  in  the  layer,  i.c.  OVC,  SCT, 
BKN,  CLR 
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TABLE  2.16  (Continued) 


hc  - Ceiling  layer  designator,  blank  if  not  a ceiling  layer 
M if  a measured  ceiling  layer 
E if  an  estimated  ceiling  layer 
MV  or  EV  if  the  ceiling  layer  is  variable 


TABLE  2.17  Format  for  Fast  Entry  of  Service  C Synoptic  Data 


Line 


1 

MMdd,  HHmm 

“ 

2 

HI,  Ty#  Nc , Nw,  VV 

3 

^ » l^h  * ^'L>  CM,  Cj-j , W 

4 

IV  W 

5 

PPP,  TT,  TdTd 

MMdd  - Four  digit  date  code,  as  described  in  Table  2.16 
HHmm  - Four  digit  time  code,  as  described  in  Table  2.16 
III  - Three  letter  station  identifier 
Ty  = C for  synoptic  data 

Nc  = Number  of  layered  cloud  groups.  This  number  was  0 tor  ail  ot 
the  synoptic  messages  in  our  data  set 

Nw  = Number  of  weathers,  equal  to  1 for  all  of  the  synoptic  menaces 
in  our  data  set 

Pl’P  - Pressure,  mb 

TT  - Surface  temperature,  °C 

TdTd  = Dew  point  temperature,  °C 

N,  Nh,  CL,  h,  CM,  CH,  VV  and  ww  are  described  in  Table  2.t> 
Consecutive  commas  are  used  to  indicate  missing  data  items. 
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TABLE  2.18  Format  for  Fast  Entry  of  Radiosonde  Data 


Data  List 


2 to  Lt  + 1 


MMdd,  HHnim,  III,  Tv 


PPPP,  hhh , TT,  TdTd 


MMdd  - Four  digit  date  code,  as  described  in  Table  2.16 

HHmm  - Four  digit  time  code,  as  described  in  Table  2.16 

III  - Three  letter  station  identifier 

Ty  = R for  Radiosonde  Data 

1^.  = Number  of  levels  of  data 

PPPP  = Level  pressure,  mb 

hhh  = Level  height,  m 

TT  = Level  temperature,  °C 

TdTd  = Level  dew  point,  °C 


Consecutive  commas  are  used  to  indicate  missing  data  items. 
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3 UTILITY  PROGRAM 


3.1  PURPOSE  AND  DESCRIPTION 

Numerous  special  purpose  or  utility  programs  .voie  t~  i ;Ci 
use  in  performing  specific  tasks  in  the  sensitivity  study  and  demonstra- 
tion run.  The  purpose  and  description  of  each  of  these  programs  togethei 
with  a listing  is  given  in  the  subsections  which  follow.  Each  subsection 
title  states  the  name  of  the  main  program  which  is  described  tneiein 
together  with  associated  subprograms.  The  system  on  which  the  piogiams 
were  implemented  is  indicated.  The  program  listings  are  given  in  the  appendix. 


3.1.1 


Program  COLUMN  (Decsystem  10) 


Program  COLUMN  reformats  the  printout  of  data  via  the  fast 
entry  format.  The  reformatting  done  by  COLUMN  results  in  the  data 
being  output  five  columns  to  a page  rather  than  one  column  as  on  input. 
This  intermediate  reformatting  was  done  to  facilitate  error  checking. 


Program  CONVE R derives  the  relative  UTM  coordinates  of  each 


station  in  the  data  sample  and  generates  the  FORTRAN  codes  used  in 
subroutine  CONVRT  which  convert  the  three  letter  station  code  to 
relative  UTM  coordinates.  Program  CONVER  utilizes  subroutine  UTM  and 
requires  the  specification  of  the  longitude  and  latitude  of  the  southwest 
corner  of  the  CFAS  window  and  the  central  meridian  of  the  window. 
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3.1.3  Program  SYNOP  (Decsystem  10) 


Program  SYNOP  converts  Service  C Synoptic  data  from  the  fast 
data  entry  format  to  the  format  readable  by  the  CFAS.  Subroutine 
CONVRT  is  called  by  S\N!OP  to  supply  the  relative  UTM  coordinates  to. 
the  stations  which  are  identified  in  the  fast  data  entry  format  by  then- 
call  letters. 


m 


- 

3.1.4  Program  AIRWAY  (Decsystem  10) 

Program  AIRWAY  converts  Service  A Airways  reports  from  the  fast 
data  entry  format  to  the  format  readable  by  the  CFAS.  Subroutines 
CLOUDS  and  WHETR  are  called  by  AIRWAY  to  convert  the  cloud  and 
weather  data  from  the  fast  data  entry  to  the  CFAS  format.  Subroutine 
CONVRT  is  called  by  AIRWAY  to  supply  the  relative  UTM  coordinates  for 

J i 

the  stations  which  are  identified  in  the  fast  data  entry  format  by  their 
call  letters. 

x 

N 
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3.1.5  Program  . NEWUTM  (Univac  1106) 


1 


1 


Program  .NEWUTM  derives  relative  UTM  coordinates  for  all 
stations  in  the  data  base.  The  latitude  and  longitude  of  the  southwest 
corner  of  the  window,  as  well  as  the  central  meridian  are  contained  in 
a data  statement.  Subroutine  .UTM  is  called  by  .NEWUTM  to  calculate 
the  UTM  coordinates.  Program  .NEWUTM  was  used  in  finding  the 
optimal  location  of  the  CFAS  window,  i.e,  a location  which  would 
maximize  the  number  of  stations  included  within  the  window  and  its 
border. 


i 
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3.1.6  Program  .CFMAIN/STAT1  (Univac  1106) 

Main  program  .CFMAIN/STAT1  and  subroutine  . EXEC1/STAT1 
which  are  modifications  of  .CFMAIN  and  EXEC1,  respectively,  were  used 
to  interpret  at  one  time  the  airways,  synop,  and  RAOB  observations  con- 
stituting the  complete  data  set  for  the  statistical  analyses  and  sensitivity 
evaluations  of  the  CFAS.  The  interpreted  observations  are  ouput  to  a 
mass  storage  file  to  be  subsequently  accessed  by  the  routines  which 
create  the  reduced  or  fractional  density  data  sets. 


3.1.7  Program  .INTPRT  (Unlvac  1106) 


Main  program  .INTPRT  reads  and  prints  out  in  a concise  and 
readable  format  a selected  mass  storage  file  of  interpreted  observations 
created  by  . CFMAIN/STAT1 , .OBSTIM,  or  .STDENS.  The  purpose  of 
.INTPRT  was  to  verify  the  contents  of  the  fractional  density  data  sets. 
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Main  program  .OBSTIM  reads  a selected  file  of  interpreted 
observations  and  sends  to  an  assigned  file  those  observations  taken 
between  specified  times.  .OBSTIM  was  used  to  create  data  sets  of 
observation  spanning  particular  time  periods. 


3.1.9  Program  . STDENS  (Univac  1106) 

Main  program  . STDENS  creates  the  fractional  density  and  silent 
area  data  sets  from  the  full  density  set.  The  procedure  employed  in 
.STDENS  is  specific  to  stations  within  the  "Alabama  Square  . The  3/4, 
1/2,  1/4,  1/8  and  1/16  full  density  and  silent  area  data  sets  are 
created  by  deleting  from  the  full  data  set  the  observations  from  specific 


stations . 


3.1.10  Program  , CFMAIN/1JPL  (Univac  11  Ob) 

Main  program  .CrMAlN/ljrL  and  subprogram  . LXCC1/IJPL  which 
are  modifications  of  .CPMAIN  and  . CXEC1,  respectively,  were  designed 
to  create  the  1FILL,  Jl'ILi;  and  BASF  file  for  each  of  the  fractional,  as 


well  as  the  full  and  silent  area  data  sets. 


3.1.11  Program  . GFMAIN/TSK3RI  (Univac  1106) 

Main  program  . CFMAIN/TSK3RI  and  subprograms  . EXEC1/TSK3RI 
and  . EXEC2/TSK3RI  are  modifications  of  . CFMAIN , .EXEC1  and  . EXEC 2 . 
The  modified  routines  were  used  to  perform  only  the  functions  involved 
in  analyzing  the  observations  and  creating  the  cloud-fog  data  base. 

This  was  done  to  facilitate  the  rapid  and  economical  evaluations  of  the 
results  obtained  with  the  fractional  and  silent  area  data  sets. 
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3.1.12  Program  .iTSTAT  (Univac  1 1 Ob) 

Main  program  . CFSTAT  roads  in  a user  selected  cloud-tog  data 

base  (CTDB)  trom  mas:  storage  and  calculates  various  statistics  of  the 

sky  cover,  coiling,  visibility,  piesent  wiathti  anil  layered  cloud  cover. 

nto  following  subprogram.,  are  called  by  . CTSTAT: 

.rrii.lNG  . S'fATPK 

.CiA'Ch  .VIMUL 

.LAYERS  . W1IATH  R 

Subprograms  .iTH.lNib,  .COVER,  .LAYERS,  .V1S11U,  and  . WEATHk 
scan,  rospoctivoly , tlio  values  ot  ceiling,  total  sky  cover,  layered  clou 
cover,  visibility  and  present  weather  stored  in  the  CFDR  to  determine 
the  numbot  of  missing  entries  and  sends  to  subprogram  . S TAT  IK  the  non 
missing  entiles  trom  which  the  statistics  are  calculated. 

Subprogram  .STATER  calculates  the  mean  value  and  variance  ot 
a specific  parameter  (i.e.  ceiling,  visibility,  etc.)  for  a selected  and 
ground  truth  CEDB.  The  covariance,  correlation  coetticnnt  and  loot 
mean  square  error  of  the  parameter  with  respect  to  ground  truth  CEDE 


are  also  calculated  by  . STATPK. 


4 


CP  AS  SENSITIVITY  ANALYSIS 


4.1  ALL  DATA  ANALYSES 

•l . I . 1 PEAS  Control  Parameter  Ettects 

•An  analysis  was  performed  tv'  determine  sensitivity  or  the  OEAS 
analysis  output  to  changes  in  usoi  inputted  oontiol  parameters.  All  data 
tor  7 March  1977  0 0 J i- 1'.  and  the  previous  1J  hours  wore  used  m the 
control  parameter  sensitivity  analysis.  Table  l . I contains  a matrix  ot 
the  control  parameters  and  values  used  for  each  ot  IS  separate  computet 
tuns.  Ihe  mound  truth  values  ot  the  control  parameters  are  identified 
separately  in  the  table.  Each  run  was  made  by  systematically  varying 
one  control  parameter  at  a time. 

Rather  than  attempt  a subjective  comparison  ot  one  map  output 
with  another  tor  the  many  runs,  a set  ot  objective  methods  was  used  tv' 
provide  a quantitative  measure  ot  "qoodness  ot  tit".  Three  such  nnlhois 
were  used,  namely,  correlation  coefficient,  root  mean  square  error  (RMSE), 
and  percentavie  of  arid  joints  missing.  Each  ot  these  methods  was  applied 
to  tin'  400  qrid  points  tor  the  analysis  window  eaeompassiiiii  Alabama. 

Only  those  qrid  points  where  an  analysis  was  possible  were  included  m 
computing  the  correlation  aiivl  RMSE.  Tlie  virnl  point  where  an  analysis 
was  not  possible  (missinq)  were  tabulated  separately  and  a pereentaqe 


of  missing  grid  points  was  computed.  Thus,  a single  value  for  each 
method  represents  a “goodness  of  fit"  for  the  entire  analyzed  map. 

Let  us  now  consider  the  control  parameter  values  used  in  the 
analyses,  as  shown  in  Table  4.1.  The  TIME  and  TYMOLD  are  constant 
values  for  the  entire  run  series  because  the  time  to  start  accepting 
observations  was  always  30  minutes  into  the  day  of  7 March  1977  (to 
include  specials  reported  at  0026Z)  and  extended  back  720  minutes  in 
time  (12  hours)  to  include  old  observations.  Two  different  patterns  were 
used  to  search  squares  for  observations  surrounding  a grid  point.  The 
ground  truth  values  had  four  possible  search  squares  to  implement,  1, 

2,  3,  or  8 grid  squares.  All  observations  within  one  grid  square  would 
first  be  used  in  deriving  an  analysis  at  the  grid  point.  If  no  observa- 
tions were  obtained,  rhe  program  would  index  to  the  second  search  square 
value  and  search  for  available  observations  to  perform  an  analysis.  This 
procedure  would  be  repeated  until  the  maximum  number  of  search  squares 
allowable  had  been  attained,  i.e.  8 for  the  ground  truth  run.  A separate 
run  was  made  using  only  two  possible  search  square  values,  i.e.  1 and 
4.  The  performance  of  this  latter  run  compared  with  ground  truth  is 
shown  in  the  bottom  row  of  Table  4.2  for  four  meteorological  parameters, 
total  sky  cover,  ceiling,  visibility,  and  present  weather,  and  for  the 
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TABLE  4.2  Results  Using  Variable  CFAS  Control  Parameters 
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percent  cloud  cover  within  tin'  nine  layers  specified  in  C FAS  to  cover 
tlie  altitude  range  from  the  surface  to  JOuO  meters.  The  correlation 
coefficient,  percent  of  missing  grid  points,  and  RMSE  values  are  pte- 
sented  in  that  order  for  each  meteorological  parameter  associated  with 
NS8Q-2,  ISSQ(lHl  and  ISSQ(2)  1.  For  sky  cover,  the  correlation  was 
0.34,  tiie  percent  missing  grid  point  data  was  2.7  and  the  RMSF  was 
8 percent.  The  units  of  the  RMSE  for  each  meteorological  variable  is 
the  same  as  those  listed  in  Table  1.1. 

Notice  that  of  all  the  runs,  the  only  time  it  was  not  possible 
to  obtain  an  analysis  at  a grid  point  was  when  only  two  search  squares  wei 
used.  A degradation  is  also  shown  via  the  correlation  coefficient  and 
RMSE  tor  all  meteorological  variables.  The  average  Univac  11  Ob  tin  to 
complete  an  analysis  was  3 minutes  for  all  data  and  -1  search  squares. 

By  having  only  two  search  values  the  computer  time  necessary  for  con  ■ 
ducting  an  analysis  was  reduced  by  18  percent.  Thus,  it  is  possible 
to  reduce  the  computer  time  significantly  by  varying  the  number  of  starch 
squares.  The  results  do  show,  however,  that  care  must  be  exercised  in 
the  cnoice(s)  of  search  square  values. 

The  effect  of  changing  the  distance  parameter,  DSP,  from  50 
to  150  km  had  the  most  pronounced  influence  on  changing  the  analysis 


output.  The  purpose  of  this  parameter  is  to  form  a single  "best  report 
on  weather  observed  by  two  or  more  very  close  stations,  those  within 


a DSP  distance  of  one  another.  The  advantage  of  this  parameter  in  the 
objective  analysis  scheme  is  that  it  allows  an  operationally  more 
important  observation  to  supersede  a less  critical  nearby  observation. 

This  procedure  is  applied  on  a variable  by  variable  nature  so  that  the 
"best  report"  could  contain  a ceiling  observation  from  one  station  and  a 
visibility  observation  from  another  close  station..  The  intent  was  to  use 
this  procedure  only  to  handle  simultaneous  observations  in  time  that 
were  also  very  close  in  space,  say  within,  plus  or  minus,  one  grid  unit 
of  the  final  analysis  mesh.  From  the  results  in  Table  4.2  it  can  be 
seen  that  increasing  DSP  to  150  km  or,  plus  or  minus,  three  grid  units 
is  too  coarse.  Although  this  parameter  has  no  effect  on  missing  grid 
points , it  does  provide  the  worst  correlation  and  largest  RMSE  for  all 
thirteen  listed  weather  variables.  It  is  clear  that  this  parameter  is 
vitally  important  to  the  quality  of  the  analysis. 

The  standard  DSP  value  of  50  proved  to  be  perfectly  acceptable 
for  the  data  dense  region  of  southeastern  U.S.A.  There  were  no  "dual" 
stations  reporting  within  50  km  of  one  another  for  this  region  so  a DSP 
of  50  was  the  best  value  to  select  to  provide  spatial  resolution  for 
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ground  truth  analyses.  A question  that  does  not  arise  with  this  normal 


type  data  sample  but  that  could  be  important  under  close  and  conflicting 
battlefield  observations  is,  what  is  the  optimum  DSP  value  that  provides 
discrimination  among  reports,  maintains  spatial  resolution,  yet  empnasi^c 
operationally  critical  features  ? 

Distance  and  time  scale  factors,  used  to  weigh  the  value  or 
influence  of  an  observation  on  a grid  point,  were  varied  to  allow  for  a 
range  of  values  from  half  as  small  to  twice  as  large  as  the  ground  truth 
value,  as  shown  in  Table  4.1.  The  first,  second,  and  tnird  aiitauce 
and  time  constants  (DIST  ( ) and  TYMC  ( ) ) appear  in  the  analysis 
scheme  to  treat  convective  (local),  convective  with  middle  clouds 
(mesoscale),  and  all  other  weather  conditions  separately.  In  general, 
the  sphere  of  influence  of  local  weather  is  limited  to  small  distances 
and  time,  i.e.  DIST  (1)  = 10,  20,  or  40  km  and  TYMC  (I)  - 2 5,  50,  or 
100  minutes.  The  range  of  the  mesoscale  and  other  distance  unu  tin,_ 
constants  is  given  in  Table  4.1.  It  can  be  seen  from  Table  4.2  that 
variations  in  the  local  and  mesoscale  cislance  and  time  constants  hove 
no  effect  upon  the  analysis  results.  The  third  distance  and  time  con- 
stants that  applied  to  all  non-convective  weather  situations  did  produce 
variations  in  the  computer  output  analyses.  However,  these  differences 


I were  extremely  small  for  all  variables  except  visibility,  which  should  be 

looked  into  further  in  view  of  the  importance  of  this  parameter  to  many 
Army  operations.  The  main  point  to  be  made  is  that  the  analysis  results 
are  not  sensitive,  nor  critical,  to  distance  and  time  weighting  factors. 

Overall  it  can  be  stated  that  no  violent  changes  in  the  analysis 
due  to  changing  input  parameters  were  detected  that  would  indicate  any 
instabilities  in  the  analysis  procedures. 
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4.2 


VARIABLE  DATA  DENSITY  ANALYSES 


4.2.1  Correlations  and  RMS  Errors 

The  ground  truth  CFAS  control  parameters  described  in  section 

4.1.1  were  applied  separately  to  every  available  station  observation 
and  then  to  variable  station  densities  to  perform  a series  of  analyses 
for  0026Z  March  7,  1977.  Maps  showing  full,  3/4,  1/2,  1/4,  1/8, 
and  1/16  station  density  distributions  are  given  in  Figures  2.12  to  2.17, 
inclusive.  The  correlation  coefficient,  RMSE,  and  percentage  of  grid 
Points  missing  an  analysis  were  derived  tor  each  station  density  using 
the  full  data  set  as  ground  truth.  These  computations  were  made  for 
each  analysis  of  the  400  inner  grid  points  that  encompass  Alat-ama.  l'hi 
correlation  coefficient  and  RMSE  computations  were  made  only  with 
respect  to  those  grid  points  where  an  analysis  was  possible.  The  per- 
centage missing  value  accounts  for  the  number  of  grid  points  where  an 
analysis  was  not  possible  even  though  the  maximum  search  square 
available  for  use  was  eight  grid  units  or  200  km.  Figure  4.1  shows 
that  the  percentage  of  missing  grid  points  remains  essentially  zero  until 
an  average  distance  of  300  km  or  more  exists  between  stations.  This 
average  distance  of  300  km  corresponds  to  only  1/8  the  number  of 
stations  where  all  are  included.  These  missing  grid  poin1  results  at" 


independent  of  the  type  of  meteoroloqic.il  variable.  It  can  also  be  seen 


from  figure  4.1,  that  the  C'l’AS  analysis  program,  using  a ma> 


iximum  ot 


eiunt  grid  unit  search  squares,  can  force  an  analysis  at  more  than  '10', 
.>1  tin  grid  points  even  tliougu  onl>  i t>  Ine  einnaal  Mniioiis  ai. 
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between  st.itions  increases.  il.ose  ic.-ailts  weie  obtained  loi  tom 
meteorological  variables,  sky  cover,  visibility,  ceiling,  and  present 
weather.  All  decay  curves  are  reasonably  well  behaved  except  tor  one 
or  two  |>oints  on  the  ceiling  display.  In  general,  as  the  avetage  dis- 
tance between  stations  doubled  from  100  to  200  km,  the  correlation 
coefficient  for  sky  cover  and  visibility  dropped  to  0.7  which  means  that 
50/.  of  the  variance  was  unexplained.  Conditions  are  oven  worse  tor 
analysing  ceiling  and  present  weather  variables  ft  out  limited  weatlui 
observations.  In  this  cast-  a doubling  ot  distance  between  stations 
resulted  in  only  about  lt>%  of  the  variability  being  explained. 

Root  mean  square  error  (KMSi;)  lesults  .ire  piesente.l  in  1’igures 
4.«  to  4.4  for  sky  cover,  visibility,  ceiling,  and  present  weather.  All 
ot  these  meteorological  variables  show  a rapid  deterioration  in  the 
KMSi:  as  distance  between  observing  stations  increases  from  100  to 
200  km.  Beyond  that  point,  deterioration  either  slows  down,  levels  oft, 


or  improves.  The  worst  RMST  is  20  percent  for  sky  cover,  1 km  for 


ceiling,  13  km  for  visibility,  and  30  category  units  for  present  weather. 
The  impact  of  data  density  on  the  quality  of  the  analysis  is  shown 
graphically  by  the  rapid  decay  or  deterioration  as  observational  densities 
decrease.  Another  factor  that  is  adversely  affected  as  station  density 
decreases  in  the  increased  running  time  required  to  maki  a computei 
analysis.  It  was  found  that  the  computer  execution  time  to  perform  an 
analysis  doubled  as  the  number  of  observation  stations  decreased  by  a 
half.  To  state  this  another  way,  less  data  icquius  longer  ceinputi  i 
runs.  The  reason  for  this  is  the  fact  that  the  computer  takes  more 
time  to  search  the  area.  This  increases  as  the  square  ot  tiu  distance 
to  the  nearest  station,  and,  in  the  process,  a fargei  number  ot  distant 
observations  are  found  that  must  be  interpreted  and  weighted  properly 
to  complete  the  analysis  at  each  grid  point.  It  is  possible  \o  u imm:.  . 
th  is  time  effect  for  any  given  data  distribution  by  aeliv  i nig  less  si.  i • ■ 
Cl’AtS  control  parameters,  such  as  the  number  and  interval  ot  search 
squares  allowable. 

Although  all  data  density  analyses,  including  silent  men,  used 
only  routine  surface  and  upper  air  observations , it  should  be  mentioned 
that  the  CPAS  has  the  built-in  capability  to  incorporate  pilot  repotts  and 
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meteorological  satellite  information.  Such  observations  could  greatly 
improve  the  available  data  density.  At  the  present  state-of-the-art. 
satellite  data  would  definitely  be  helpful  in  improving  sky  cover  and. 
to  a far  lesser  extent,  present  weather,  ceiling,  and  visibility  result 
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Figure  4.11 


Silent  Area  Grid  Point  Krror  Values 
For  Visibility  (Kilometers) 


grid  point  analyses  in  the  southeast  portion  of  Figures  4.10  to  4.12. 

The  maximum  number  of  search  squares  allowable  for  this  study 
was  eight.  Increasing  the  number  of  search  squares  to  ten  for  the 
northeast  sector  and  to  thirteen  for  the  southeast  sector  of  the  map 
would  have  forced  an  analysis  at  all  grid  points  so  that  there  would  be 
no  missing  analysis  points.  The  CFAS  search  control  parameter  can 
always  be  set  as  the  greatest  distance  between  a grid  point  and  the 
nearest  observation  to  insure  no  missing  grid  point  analyses  even  in 
silent  areas.  The  further  a grid  point  is  removed  from  a data  point  and 
the  smaller  the  scale  of  the  phenomena  the  more  questionable  the 
analysis.  On  the  other  hand,  some  have  argued  that  such  an  analysis 
expanded  to  include  climatological  and  forecast  or  extrapolated  informa- 
tion would  produce  more  reliable  estimates  than  a pilot's  guess  of 
silent  area  conditions.  The  CFAS  can  accept  forecast  information  as 
data  inputs  to  its  objective  analysis.  The  CFAS  does  not  presently 
have,  however,  the  capability  to  perform  its  own  trend  or  forecast 
analysis  for  extrapolating  information  into  a silent  area. 

It  can  be  seen  from  Figure  4.10  that  the  CFAS  sky  cover 
analysis  not  only  has  no  errors  in  the  data  region  but  also  does  re- 
markably well  in  the  silent  area.  It  is  interesting  to  note  that  neither 
statement  can  be  made  concerning  analyses  of  mesoscale  features 


123 


associated  with  visibility  or  base  of  the  lowest  cloud,  shown  in  4.12 
and  4.11,  respectively.  In  both  of  the  latter  cases,  lack  of  stations 
in  the  southern  part  of  the  data  area  adversely  effects  analyses  in  both 
the  silent  and  data  area.  It  can  be  stated  that,  in  general,  for  all 
three  weather  variables,  results  deteriorate  as  distance  increases  into 
the  silent  area.  In  fact,  the  largest  visibility  error  exceeds  1 0 km  an  I 
the  largest  error  in  the  lowest  cloud  base  exceeds  3 km.  The  importation 
of  having  some  observations,  especially  of  mesoscale  features,  can  not 
be  overemphasized,  and,  as  discussed  previously,  reconnaissance  or 
satellite  reports  would  be  valuable  additions  to  silent  area  analyses. 

This  is  further  emphasized  by  the  "goodness  of  fit"  results  shown  in 
Table  4.3.  In  this  table  the  correlation  coefficient,  percent  of  missing 
grid  point  analyses,  and  RMSE  are  tabulated  for  nine  meteorological 
variable  analyses  completed  within  the  500  km  square  region  that  in- 
cludes half  a data  rich  and  half  a silent  area.  A maximum  of  18.5* 
of  the  grid  points  were  missing  an  analysis  for  the  first  cloud  layer, 
whereas  all  other  variables  had  a 16.5%  missing  data  analysis  rate. 

The  correlation  coefficients  and  RMSE  values  were  the  lowest  and 
highest,  respectively,  for  such  weather  variables  as  ceiling  and  layered 
cloud  cover.  A correlation  coefficient  of  only  0.31  was  obtained  for 
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TABLE  4.3  Silent  Area  Goodness  of  Fit  Values 


Correlation 

Coefficient 


% 

Missing 
Grid  Points 


RMSE 


p 


the  percent  cloud  cover  in  the  first  layer  which  extends  from  the  surface 
to  45  meters  above  ground  level  (AGL).  The  correlation  Improved  for 
percent  cloud  cover  within  the  second  layer  (from  45  to  91  meters  AGL) 
and  essentially  continued  to  improve  for  each  successively  higher  layer. 
Thus,  the  greatest  difficulty  in  analyzing  percent  cloud  cover  within 
particular  layers  occurred  in  layers  closest  to  the  ground.  It  is  pre- 
cisely in  this  region  that  there  were  fragmented  clouds  and  low  scud  which 
drift  rapidly  into  and  out  of  sight  of  the  observer  and  present  the 
worst  objective  analysis  results. 
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5 ARMY  AVIATION  WEATHER 
DEMONSTRATION  MATERIAL 


5.1  MAP  DISCUSSION  OF  RESULTS 

The  purpose  of  this  section  is  to  demonstrate  a variety  of 
capabilities  that  are  possible  using  computer  applications  to  Army 
aviation  weather  problems.  In  order  to  appreciate  some  of  the  problems, 
a typical  teletype  weather  report  is  reproduced  in  Fig.  5.1.  The  problem 
that  a pilot  encounters  in  attempting  to  acquaint  himself  with  the 
current  weather  is  compounded  by  the  fact  that  weather  information  is 
coded  primarily  for  meteorologists,  requires  a knowledge  of  local 
station  call  letters  and/or  numbers,  is  filed  on  clip  boards  according  to 
FAA  transmission  circuit  number,  or  appears  on  difficult  to  read  facsimile 
weather  charts.  In  order  for  a pilot  to  obtain  the  latest  weather  signifi- 
cant to  his  flight  it  is  necessary  for  him  to  screen  or  sift  through  all 
special  weather  observations  to  obtain  and  decifer  those  important  to 
his  mission.  Often  on  urgent  missions  and  in  the  hustle  of  immediate 
operations  it  is  not  possible  for  a pilot  to  maintain  proficiency  in  his 
primary  function  while  maintaining  cognizance  of  weather  restrictions  on 
his  performance. 

The  present  Army  weather  support  system  evolved  from  ex- 
tension of  the  labor  intensive  concept  of  having  trained  meteorologists 
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Fig.  5.1  Typically  Difficult  to  Read  Teletype  Weather  Report 
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providing  pilot  weather  briefings.  The  time  period  of  the  1980's  and  the 
expected  available  trained  personnel  and  support  logistics  suggest  that 
an  automated  method  of  providing  weather  information  can  not  only 
provide  better,  up  to  date,  and  operationally  tailored  information,  but 
can  do  so  in  a manner  that  greatly  reduces  manpower  while  increasing 
efficiency  within  a battlefield  environment.' 

A pilot  controlled  and  operated  computer  weather  system  can 
provide  upon  his  command  not  only  the  latest  update  of  weather  informa- 
tion tailored  to  his  needs  but  can  also  provide  the  past  sequence  of 
events  and,  if  necessary,  can  provide  a short-time  projection  of  future 
weather  events . 

For  any  geographical  or  weather  analysis  region,  the  computer 
can  provide  such  operationally  important  features  as  the  topography  and 
geopolitical  boundaries  upon  which  an  analysis  can  be  made.  For  our 
purposes.  Fig.  5.2  shows  these  features  for  the  Alabama  region  where 
we  will  focus  our  attention  for  this  study  to  demonstrate  some  of  the 
products  that  are  possible  to  support  Army  operations. 

Before  proceeding  to  discuss  the  many  figures  that  describe 
the  computers  capability  to  present  time  series  of  weather  events  or 
many  weather  events  for  one  time  period,  it  should  be  mentioned  that 
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this  is  i very  limited  sample  of  what  is  possible.  It  should  also  be 
emphasized  that  in  order  for  the  maximum  benefit  to  be  achieved  from 
automated  weather  applications  to  Army  operations,  it  is  necessary  that 
both  the  Army  operational  and  meteorological  community  identify  their 
mutual  problems,  assign  priorities,  and  work  collectively  to  attain 
defined  goals. 

We  have  selected  a time  series  of  meteorological  events 
that  would  be  of  interest  to  aviation  personnel.  Computer  analyses 
were  made  at  odd  times  to  be  more  characteristic  of  pilots  interrogating 
the  computer  terminal  on  a random  basis.  The  time  series  analyses  are 
given  for  five  meteorological  variables.  Figs.  5.3  to  5.7  present  com- 
puter analyses  of  per  cent  sky  cover  for  five  separate  times  beginning 
2130Z  26  February  1977  and  ending  0640Z  27  February  1977.  Figs.  5.8 
to  5.11  present  the  height  of  the  lowest  cloud  base  beginning  G040Z 
and  ending  0640Z  27  February  1977  . Figs.  5.12  to  5.15  present  the 
ceiling  height  and  Figs.  5.16  to  5.19  depict  the  horizontal  visibility 
at  the  surface  for  the  same  time  interval  as  the  lowest  cloud  base 
analyses.  Fig.  5.20  to  5.24  also  cover  the  same  time  interval,  plus 
an  additional  analysis  of  present  weather  conditions  at  07  20Z 
27  February  1977.  A color  code  was  used  for  all  the  analyses  to  depict 
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categories  of  weather  conditions.  In  all  cases,  red  was  used  to  indicate 
the  worst  aviation  weather  conditions,  then  came  yellow,  green,  and 
blue,  in  that  order,  with  white  areas  representing  the  best  flight  weather. 

Sky  cover  analysis  for  2130Z  26  February  in  Fig.  5.3  shows 
the  series  begins  with  most  of  the  area  having  clear  or  partly  cloudy 
skies.  Only  in  the  southwest  and  northwest  corners  do  broken  clouds 
exist.  Four  hours  later.  Fig.  5.4,  spotted  regions  on  the  map  have 
overcast  conditions  but  as  time  progresses,  the  entire  area,  as  seen  in 
Fig.  5.7,  becomes  engulfed  by  nearly  an  unbroken  solid  mass  of  clouds. 

The  lowest  cloud  base  analysis  series  shows  the  detail 
character  that  is  attainable  with  the  computer  system.  This  is  one 
parameter  that  is  difficult  to  obtain  a mental  picture  of  what  exists  by 
reading  teletype  reports.  It  is  also  a parameter  that  is  not  routinely 
analyzed  and  transmitted  on  the  facsimile  circuit.  Yet  it  is  a parameter 
that  has  value,  especially  to  helicopter  pilots  that  do  not  have 
sophisticated  instruments  to  fly  in  clouds.  Notice  in  Fig.  5.8,  there 
are  two  small  regions  where  the  lowest  cloud  bases  are  below  0.15  km. 
The  corresponding  sky  cover  map  , Fig.  5.4,  shows  overcast  skies  for 
these  same  areas  but  the  ceiling  map.  Fig.  5.12,  shows  the  overcast 
skies  are  above  1.2  km.  A corresponding  map.  Fig.  5.16,  shows 
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visibility  exceeds  4 km  for  the  same  geographical  location  as  those  low 
clouds  while  the  present  weather  map.  Fig,  5.20,  shows  fair  for  one 
and  rain  for  the  other  location.  Thus,  by  looking  at  several  of  these 
maps  for  the  same  time  period  one  can  obtain  a good  mental  picture  of 
weather  conditions  prevailing  at  a given  location.  This  feature  will  be 
demonstrated  more  in  the  later  discussion  on  Figs.  5.25  to  5.33.  As 
far  as  the  lowest  cloud  base  map  series  is  concerned,  the  patchy 
nature  and  detail  of  this  variable  is  clearly  illustrated  along  with  its 
rapid  time  changes  as  shown  in  Fig.  5.8  to  5.11.  Field  Army  flights 
are  frequent  but  usually  less  than  an  hour  in  duration.  Having  a 
system  that  keeps  a pilot  "weatherwise"  of  the  latest  details  would 
prove  valuable  in  planning  and  executing  a mission. 

Categories  desired  for  displaying  an  analyzed  variable  can 
be  inputted  by  the  user  to  provide  flexibility  in  satisfying  requirements 
of  different  missions.  Since  most  Army  flights  are  flown  at  altitudes 
below  1.2  km,  the  ceiling  and  lowest  cloud  height  categories  were 
chosen  to  provide  analysis  resolution  in  the  lowest  layers.  A cloud 
ceiling  is  only  reported  for  those  clouds  that  cover  a sufficient  portion 
of  the  sky  to  constitute  either  broken  or  overcast  conditions.  For  this 
reason  considerable  differences  exist  between  the  lowest  cloud  base 
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maps  and  the  ceiling  time  series  maps  in  Figs.  5.12  to  5.15.  Very 
low  ceilings,  below  0.15  km,  remain  stationary  in  the  southern  part  of 
Alabama  while  very  low  ceilings  move  systematically  from  west  to  east 
across  the  northern  part  of  the  state.  Horizontal  visibility  begins  the 
series  exceeding  A km  for  all  but  a small  (25  to  50  km)  portion  of  the 
map,  Fig.  5.16.  Only  the  southern  part  of  the  state  experiences 
visibilities  less  than  l km  and  for  two  different  reasons.  A small  strip 
of  low  visibilities  produced  by  fog  remains  static  in  the  southeast 
corner  of  the  state.  Those  low  visibilities  in  the  southwestern  corner 
of  the  map,  Fig.  5.19,  are  produced  by  rain  from  thunderstorms,  as 
shown  in  Fig.  5.23.  At  a glance  it  is  possible  to  delineate  regions  of 
high  visibility  and  to  track  the  movement  or  stagnation  of  visibilities 
within  operationally  defined  categories. 

Present  weather  analyses  begin  with  a touch  ot  snow,  rain, 
and  showers  or  thunderstorms.  Fig.  5.20,  but  most  of  Alabama  is  under 
the  influence  of  fair  weather.  Three  hours  later,  in  Fig.  5.21,  the 
snow  in  the  north  has  changed  to  rain,  the  rain  in  the  south  has  stopped, 
thunderstorms  remain  in  the  west  but  begin  to  appear  in  the  north  ahead 
of  the  steady  rain,  and  radiational  cooling  under  fair  skies  has  resulted 
in  fog  conditions  in  the  southeast  corner.  An  hour  and  half  later. 

Fig.  5.22,  the  radiation  fog  remains,  the  rain  returns  in  the  south  and 
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expands  in  the  north,  and  showers  and  thunderstorms  ring  the  northern 
rain  area.  After  another  hour  and  a half.  Fig.  5.23,  rain  and  showers 
progress  southeastward  and  the  radiation  fog  remains  entrenched  in  the 
southeast  corner.  Only  40  minutes  later,  Fig.  5.24,  a 200  km  wide 
band  of  thunderstorms  exists  about  a northeast-southwest  axis  extending 
from  C11A  (Chattanooga , TN)  to  MOB  (Mobile,  Al.)  while  the  southeast 
radiation  fog  continues  to  persist.  Computer  analyses  of  the  latest 
present  weather  phenomena  should  find  great  acceptance  by  pilots  and 
prove  to  be  very  beneficial. 

At  any  one  instant  in  time,  a pilot  may  wish  to  view  all 
weather  variables  that  exist  within  his  area  of  interest.  Such  a set  is 
provided  in  Figs.  5.25  to  5.33  for  nine  weather  parameters.  The  entire 
state  of  Alabama  has  overcast  skies,  Fig.  5.25,  with  a good  portion  ot 
it  containing  low  clouds,  Fig.  5.2b,  whose  bases  are  less  than  0.3  km 
Although  all  of  Alabama  has  a cloud  ceiling,  only  the  northern  and 
southern  portions  have  ceilings  below  O.b  km,  Fig.  5.2/.  It  is 
interesting  to  note  that  only  southern  portions  experience  low  visibilitie 
Fig.  5.28,  produced  by  fog  in  the  east,  rain  in  the  midsection,  and 
showers  in  the  western  part  of  the  state,  Fig.  5.28. 

A pilot  may  wisli  to  have  information  tailored  to  his  flight 
altitude.  One  such  output  is  presented  to  show  the  extent  that  a 


particular  flight  layer  is  obscured  by  clouds.  Figs.  5.30,  5.31  , 5.32, 
and  5.33  display  the  per  cent  cloud  cover  as  overcast,  broken,  partly 
cloudy,  or  clear  for  the  flight  layers  from  surface  to  45  m,  45  to  91  m, 
91  to  183  m,  and  183  to  305  m AGL,  respectively.  Notice  the  detail 
and  structure  that  is  attainable  via  the  computer  and  how  a pilot  can 
quickly  assimilate  his  ability  to  maintain  visual  contact  within  each 
layer. 

A number  of  other  candidate  computer  output  products  include 
such  items  as: 

1 . severe  weather  tracks 

2.  areas  below  minimum  ceilings  and/or  visibilities 

3.  VFR  and  IFR  flight  paths 

4.  vertical  sections  of  flight  plans 

5.  air  to  ground  clear  lines  of  sight 

6.  best  flight  plan  for  not  being  seen  from  ground 

7.  best  flight  plan  for  photo  recon  mission 

8.  plus  many  others  where  weather  factors  influence  Army 

operations . 
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Fig.  5.3  Sky  Cover  2130Z  26  Feb  1977  Computer  Analysis 
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Fig.  5.4  Sky  Cover  0040Z  27  Feb  1977  Computer  Analysis 
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Fig*  5,5  Sky  Cover  0340Z  27  Feb  1977  Computer  Analysis 
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Fig.  5.6  Sky  Cover  0510Z  27  Feb  1977  Computer  Analysis 
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Fig.  5.8  Lowest  Cloud  Base  0040Z  27  Feb  1977  Computer  Analysi 
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Fig.  5.9  Lowest  Cloud  Base  0340Z  2.1  Feb  1977  Computer  Analysis 
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Fig.  5.11  Lowest  Cloud  Base  0640Z  27  Feb  1977  Computer  Analysis 
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Fig.  5.13  Ceiling  0340Z  27  Feb  1977  Computer  Analyst 
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Fig.  5.14  Ceiling  0510Z  27  Feb  1977  Computer  Analysis 
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Fig,  5.15  Ceiling  0640Z  27  Feb  1977  Computer  Analysis 
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Fig.  5.17  Visibility  0340Z  27  Feb  1977  Computer  Analyst 
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Fig.  5.18  Visibility  0510Z  27  Feb  1977  Computer  Analysis 
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Fig.  5.19  Visibility  0640Z  27  Tcb  1977  Computer  Analy 
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Fig.  5.20  Present  Weather  0040Z  27  Feb  1977  Computer  Analysis 


154 


white  FAIR 

SNOW 

mmm-  SHOWERS  & 

THUNDERSTORMS 


Fig.  5.22  Present  Weather  0510Z  27  Feb  1977  Computer  Analysis 
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Fig.  5.24  Present  Weather  0720Z  27  Feb  1977  Computer  Analysis 


i 

■4uhmmmmb 

Pitt 


white 


CLEAR 

PARTLY  CLOUDY 

BROKEN 

OVERCAST 


Fig.  5.25  Sky  Cover  0640Z  27  Feb  1977  Computer  Analy 
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Fig.  5.32  Cloud  Cover,  91  to  183  m AGL  Layer,  0820Z  27  Feb  197  7 
Computer  Analysis 
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Fig.  5.33  Cloud  Cover,  183  to  305  m AGL  Layer,  0820Z  27  Feb  1977 
Computer  Analysis 


6 RESULTANT  MODIFICATIONS  OF 
C FAS -C FAR  PROGRAM  ELEMENTS 


1 

I 

6.1  SUMMARY 

In  the  course  of  this  effort,  the  CFAS-CFAR  was  used  extensively 
on  a large  data  base  exhibiting  a wide  range  of  cloud  formations,  surface 
visibilities  and  weather  phenomenon.  As  a result  of  this  activity,  the 
need  for  modifications  to  several  CFAS-CFAR  program  elements  was  identi- 
fied. These  modifications  were  prompted  by  one  or  more  of  the  following 
reasons: 

1)  Results  of  the  sensitivity  analysis 

2)  Discovery  of  errors  or  bugs 

j 

3)  Improved  operation  of  the  system 

| 

The  program  elements  affected  and  the  modifications  are  described  in  the 
subsections  which  follow. 

6.1.1  FORTRAN  PROC  .BASE  (Univac  1106) 

In  order  to  minimize  the  data  storage  and  instruction  code  size  ] 

requirements  of  the  CFAS-CFAR,  the  dimensions  of  several  arrays  con- 
tained in  COMMON  (BASE)  must  be  tailored  to  the  number  of  stations  or 
reporting  sites  within  the  CFAS  window  and  border  and  the  frequency  of 
the  reports.  In  order  to  facilitate  the  changes  that  are  required  when 


the  window  is  relocated  COMMON  (BASE)  and  an  associated  DATA  state- 


ment have  been  placed  in  FORTRAN  PROC  .BASE.  This  permits  the 


incorporation  of  these  declaratory  statements  contained  in  .BASE  into  the 
several  subprograms  which  reference  them  through  the  use  of  a simple 


INCLUDE  statement.  The  affected  subprograms  are: 


. BEGIN 
. EXEC1 
.FIN  DIB 
.GETOB1 
.GET1BW 


.GET1FW 
. ITOJ 
.NOSECT 
.STOREC 
.SECTOR 


6.1.2  Main  Program  .CFMAIN/DSK  (Univac  1106) 

Main  program  .CFMAIN/DSK  which  runs  the  CFAS  was  revised 
to  incorporate  the  findings  of  sensitivity  analysis.  The  analysis  control 
parameters  whose  values  have  now  been  fixed  as  a result  of  this  study 
are  now  contained  within  DATA  statements  and  are  no  longer  input  via 
READ  statements.  The  revised  input  data  stream  now  required  for  TASK  3 
is  given  in  Table  6.1.  In  addition  to  these  changes,  there  has  also 
been  incorporated  into  .CFMAIN/DSK  a provision  to  allow  the  interpreted 
observations  to  be  printed  out  following  a TASK  2.  This  feature  is 
actuated  by  inserting  a SETC  command  in  the  runstream  at  some  point 
prior  to  the  XQT  command.  The  printout  of  the  interpreted  observations 
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TABLE  6.1  Input  Data  Stream  for  CFAS  TASK  3 


Line 

1 

2 

3 

5 

Data  Element 
TASK 

NPRT 

NOWTYM 

NBKOUT 

TIME 

TYMOLD 

I DENT 


Data  Elements  Format 

TASK,  NPRT,  NOWTYM  3110 

NBKOUT  110 

TIME,  TYMOLD  2110 

I DENT  10A6 


Definition 

TASK=3  for  an  analysis  at  all  grid  points  in  the 
window 

NPRT  0 the  grid  point  analysis  is  printed  out 
NPRT=0  no  printout  of  the  grid  point  analysis 
The  current  clock  time 

The  number  of  the  most  recent  record  on  the 
CFDB  file 

Time  of  the  most  recent  observation  to  be  used  in 
the  analysis 

Time  of  the  oldest  observation  to  be  used  in  the 
analysis 

A user  defined  sixty  character  label  which  identitu 
the  CFDB  record 


is  independent  of  the  optional  echo  printout  of  the  observations  on  read- 
in  which  has  always  existed  in  the  CFAS. 

6.1.3  Subprogram  . CFMAP  (Univac  1106) 

A minor  bug  was  discovered  in  CFMAP  relating  to  situations  in 
which  the  nearest  observation  to  a grid  point  was  nevertheless  distant 
enough  so  as  to  produce  a very  small  weighting  function.  The  modifi- 
cation incorporated  into  CFMAP  established  a minimum  value  for  the 
smallest  weighting  function.  If  the  smallest  weighting  function  is  less 
than  the  minimum  value,  a missing  value  is  inserted  at  the  grid  point 
for  the  parameter  being  analyzed. 


6.1.4  Subprogram  .CQMOBR  (Univac  1106) 

The  criteria  for  ranking  observations  in  the  building  of  a best 
report  were  modified  and  reordered.  The  criteria  for  ranking  observations 
used  in  creating  a best  report  are  now  applied  in  the  following  order: 

1 - Time  of  observation 

2 - Type  of  observation 

1 - Airways 

2 - Metar 

3 - Synop 

4 - Raob 

5 - AFG WC-3DNEPH 

3 - Value  of  observation 

4 - Urgency  of  observation 

5 - Distance  from  best  report  site 


6.1.5  Subprogram  . CXEC1  (Univac  1106) 

Modifications  of  subprogram  . EXEC1  were  made  in  connection 
with  the  use  of  PROC  .BASE  (Section  6.1.1)  and  the  optional  printout 
of  interpreted  observation  (Section  6.1.2). 


6.1.6  Subprogram  . EXEC 2 (Univac  1106) 


Subprogram  .EXEC2  was  modified  to  permit  the  incorporation  of 
the  grid  point  altitudes  via  an  INCLUDE  statement.  The  grid  point 
altitudes  are  carried  in  FORTRAN  PROC  . GRDPOIN  TALT. 

6.1.7  FORTRAN  PROC  .GRDPOINTALT  (Univac  1106) 

The  PROC  .GRDPOINTALT  contains  the  DIMENSION  and  DATA 
statements  in  which  the  array  of  grid  point  altitudes  are  stored.  As  is 
the  case  with  PROC  .BASE,  PROC  .GRDPOINTALT  is  dependent  upon  the 
location  of  the  CFAS  window.  The  use  of  PROC  elements  and  the 
INCLUDE  statements  which  are  available  on  high  level  FORTRAN'S,  such 
as  FORTRAN  V,  provide  a convenient  and  relatively  error  safe  means  of 
entering  semi-permanent  parameters  and  constants  into  the  CFAS  elements. 

6.1.8  Subprogram  .SFDINT  (Univac  1106) 

A modification  was  made  to  subprogram  .SFDINT  to  incorporate 
a test  and  rectification,  if  necessary,  to  insure  that  the  maximum  heights 
of  cloud  tops  which  are  derived  from  weather  information  are  equal  to  or 
greater  than  the  observed  minimum  bases  of  clouds. 


6.1.9  Subprogram  .SHADE 

An  error  was  uncovered  in  subprogram  .SHADE  which  resulted  in 
the  incorrect  shading  on  most  significant  present  weather  maps  of  weather 
categories  40  - 49.  The  error  was  corrected. 
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7 CONCLUSIONS  AND 
RLC  OM  M EN  D ATI  ON  S 


1.  Reduced  station  density  adversely  affects  the  accuracy  of  the 
analysis  at  individual  grid  points  since  more  distant  observations 
must  be  used  in  analyzing  the  CPAS  parameters  at  the  grid  point. 

2.  The  worst  objective  analysis  results  were  as  expected,  associated 
with  regions  having  no  data  (silent  area)  and  with  small  space  and 
time  scale  phenomena,  such  a rapidly  moving  scud  type  cloud 
near  the  earth's  surface. 

3.  The  inclusion  of  older  observations  at  stations  at  which  there 
exists  current  observations  is  of  value  only  to  the  extent  that 
missing  elements  in  the  current  observations  can  Lie  supplemented 
by  non-missing  values  in  the  older  observations. 

4.  Reduced  station  density  causes  an  increase  in  running  time.  This 
is  due  to  the  fact  that  the  CPAS  is  designed  to  systematically 
search  square  areas  of  increasing  size  around  the  grid  point  until 
a minimum  number  of  observations  is  found.  The  time  required 
for  the  search  is  inversely  proportional  to  the  density  of  the 
observations.  The  time  impact  of  observation  density  will  have  to 
be  carefully  considered  when  a computer  is  selected  for  operationally 
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implementing  CFAS.  It  is  also  possible  to  modify  and  optimize  the 
CFAS  for  specific  observational  densities  and  distributions.  This 
should  also  be  considered  before  implementation  of  the  system. 

5.  It  is  possible  to  trade  off  analysis  accuracy  against  computer  time 
required  for  an  analysis  by  varying  the  allowable  search  square 
control  parameter. 

6.  Consideration  should  be  given  to  handling  a problem  likely  to  exist 
within  a battlefield  environment  where  close  and  conflicting  observa- 
tions dictate  the  need  for  a technique  that  provides  discrimination 
among  reports,  maintains  spatial  resolution,  yet  emphasizes 
operationally  critical  features, 

7.  The  distance  and  time  factors  used  to  weight  the  influence  of 
distant  observations  to  a grid  point  value  proved  totally  adequate 
to  handle  all  meteorological  variables  tested,  with  the  possible 
exception  of  visibility  which  should  be  explored  more  thoiouyhly 
in  view  of  its  importance  to  many  Army  operations. 

8.  No  instabilities  were  detected  in  the  CFAS  analysis  procedures. 

9.  A most  significant  result  of  this  study  was  the  establishment  of  the 
fact  that  the  CFAS/CFAR  objective  analysis  and  display  techniques 
remained  stable  and  continued  to  perform  under  excessive  variations 
of  control  parameters  and  type,  distribution,  and  density  of 
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observations.  This  made  it  possible  to  fix  control  parameters  and 
eliminate  these  as  user  inputs.  The  computer  programs  have  been 
modified  to  incorporate  these  advances  and  a Univac  1100  series 
compatible  tape  containing  the  new  program  has  been  delivered  to 
Atmospheric  Sciences  Laboratory. 

10.  Presentations  of  the  CPAS  and  CFAR  analyses  should  be  made  to 
Army  personnel  to  demonstrate  some  possible  types  of  output,  to 
solicit  critical  comments  and  support,  and  to  generate  a list  of 
additional  products  that  would  increase  the  value  to  Army  users. 


PROGRAM  COLUMN 
DIMENSION  F’AGE<55,5»20) 

DOUBLE  PRECISION  INF  IL.E  , OUTF IL 
LOGICAL  EOF 
EOF '=. FALSE. 

TYPE  1 

1 FORMAT  ('  INPUT  INFILE .EXT  > 1 NUN IT » 1 OUT. EXT »OUTUNT ' ) 
ACCEPT  *,  INFILE » 1 NUN. IT  > 0UTF1L  , 1 OUT 

OPEN  < UN I T= I NUN IT »F'ILE=INFILE , ACCESS=  ' SEQIN ' ) 

OPEN  ( UN  I T=  I OUT  »FILE=QUTFIL  , ACCESS-' ' SEQQUT  ' ) 

5 DO  2 1=1 , 5 
DO  2 J=1 , 55 
DO  2 K=1 , 20 

2 PAGE( J, I ,K>='  ' 


DO  10  1 = 1,5 

DO  20  J=1 , 55 

READ  ( INUNIT »30>END=31 ) < PAGE ( J , 1 , K ) , K = 1 , 20 ) 
20  CONTINUE 

10  CONTINUE 

30  FORMAT  <20A1> 

35  DO  40  J=1 , 55 

WRITE  < IOUT » 60 ) < < PAGE < J, I »K ) , K=1 , 20 ) , 1 = 1 , 5 > 

40  CONTINUE 

WRITE  (IOUT, 11) 

11  FORMAT  ( 16X , 100 ('#')) 

60  FORMAT  ( IX , 3 ( 20A1 , 8X > , 20A1 , 7X , 20A1 > 

IF  (.NOT.  EOF)  GOTO  5 
GOTO  99 

31  EOF= . TRUE . 

GOTO  35 

99  CONTINUE 
END 


PROGRAM  CONNER 

INTEGER  ALT , EST  > NTH , LAB , DEG, DEG1 , MIN , MINI ,NEXLAB 

REAL  REFLON , REEL AT , CMRD, REFEST, REFNTH, LON. LAT, EAST r NORTH 

OPEN  (UNIT=20,FILE=' JACK.DAT' ) 

OPEN  < UN IT-22,  FILE- ' RAPP . DAT ' ) 

TYPE  1 

.1.  FORMAT  ('  REFERENCE  LONG') 

ACCEPT* , REFLON 
TYPE  2 

2 FORMA  I ( ' REFLAT ' ) 

ACCEPT* , REFLAT 
TYPE  3 

3 FORMAT  ('  CMRD ' ) 

ACCEPT* ,CMRD 

WRITE  (5,4)  REFLON, REFLAT, CMRD 

4 FORMAT  (3F10.3) 

LABEL-0 

CALL  UTM  < REFLON , REF  LAI  , REFEST , REFNTH , CMRD ) 

READ  (20,20)  STATIN, DEG, MIN, DEG1 , MINI , ALT 
LON=FL.OAT  ( DEG ) f ( FLOAT  ( MIN  ) /60  . ) 

LAT-FL0AT(DEG1 > + ( FLOAT (MINI ) /AO.  ) 

CALL  UTM  (LON, LAT, EAST, NORTH, CMRD) 

ESI- I NT ( 1000 . * ( EAST  REFEST ) ) 

NTH- 1 NT  ( 1.000  ♦ * ( NORTH -REFNTH ) ) 

WRITE  (22,30)  STATIN, EST, NTH, ALT 
10  READ  (20, 20, END-50)  STATIN, DEG, MIN, DEG1 , MINI , ALT 
L A B E L - L A B E L F 1 0 

LON-FLOAT ( DEG ) F ( FLOAT ( MIN ) /60 . ) 

I..  AT -FLOAT  ( DEG1 ) F ( FLOAT  ( MINI ) /60 • ) 

CALL  UTM  ( LON, LAT, EAST, NORTH, CMRD) 

EST ^1 NT ( 1000 . * ( EAST-REFEST ) ) 

NTH-INT ( 1000 , * ( NORTH- REFNTH ) ) 

N E X 1.  A B - L A B E L F 1 0 

WR I TE  ( 22 , 40 ) LABEL , STATIN, NEXLAB , EST , NTH , ALT 
GOTO  10 

20  FORMAT  (A, 51) 

30  FORMAT  (6X,17HIF  (STATIN  ,NE.  '^3,100')  GOTO  10,/,10X,  ' JX="  , 15, 
■f  / , 10X,  ' JY= ' , 15 ,/,  10X , ' JZ= " , 15 , / , 10X , "GOTO  500') 

40  FORMAT  (2X, 13, IX, 17HIF  (STATIN  .NE.  ',A3,7H'  GOTO  ,I3,/,10X, ' JX- 
F 15,/, 10X, ' JY=' , IS,/, 10X, ' JZ=' , 15,/, 10X, 'GOTO  500' ) 

50  CONTINUE 
END 

SUBROUTINE  UTM ( LON , LAT , EAST , NORTH , CMRD ) 

REAL  LAT, LON, NORTH 
A-63. 782064 
ARE D -6 3 . 350345 
E-. 0068 147849 
0~.017453292*LAT 
P=3600.*( CMRD -LON) 

C-COS(O) 

S=SIN(Q) 


1H0 


T=S/C 
S2-2.*S*C 
D=1 . ~<2.*S*S> 

84-2 . *S2*D 

RHO -A/SQRT ( 1 . - < 6 . 7686580E 03*S*S ) ) 

D=Q+ ( . 005076492* (Q-< . 5*S2 ) ) ) f 
f (4.29513E  05*<  < 1.5*0)  S2+(S4/8. ) ) ) 

XN1=ARED*D 

D~C*S*1 . 1752215E-1 1*P*F' 

D=D+ ( < C**3  > *S*2 . 3015189E- 23* ( P**4 ) * 
f <5.-<T*T)+<9.*<  < E*C ) **2 ) ) f (4.*< <E*C)**4) ) ) ) 
N0RTH=.9996*<XN1+<D*RH0> ) 

D=C*4 . 848 1 368E- 06*P 

D==D+<  (C**3)*(  1 .-(T*T)  f < <E+C>**2>  >*1 . 89921 15E  1 7* < P**3 ) > 
EAST= < RH0*D* . 9996 ) +5 
RETURN 
END 


r 
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PROGRAM  SYNOP 

INTEGER  INUNIT  - OUTUNT  . TESTER . DATE  r TIME  . LAYERS , UIETHRS  . VISIBL  . 

+ M ISSNG  . N . NH . CL.  . H » CM  . CH  . W . WW  . PPP  . TT  . TDTD . JX  . JY , JZ * TYPES 

DOUBLE  PRECISION  INFILE »OUTF I L 
DIMENSION  TEST ( 7 ) 


TYPE  10 

10  FORMAT  ('  NAME  THE  INPUT  UNIT  AND  INPUT  FILE. EXT') 
ACCEPT  *.  INUNIT . INFILE 

TYPE  20 

20  FORMAT  <'  NAME  THE  OUTPUT  UNIT  AND  OUTPUT  FILE. EXT  7 ) 
ACCEPT  *,  OUTUNT .OUTFIL 

OPEN  (UNI  1 - INUNIT » FILE-  INFILE  * ACCESS- 7SEQIN 7 > 
OPEN  (UNIT-  OUTUNT.  FILE=  OUTFILv  ACCESS® 'SEQOUT ' ) 

M ISSNG®  -32768 
TYPEC®  3 


30  READ  (INUNITf 100»  END®  90.  ERR®  40)  TESTER 
REREAD  110.  DATE. TIME 
READ  (INUNIT. 120)  ( TEST ( I ) . 1-1 . 5 ) 

40  REREAD  130.  STATIN.  TYPE.  LAYERS.  WETHRS.  VISIBL 
REREAD  120.  ( TEST ( I ) » 1=1 .5) 

IF  ( TEST ( 5 ) .EQ.  ' ')  VISIBL ®M I SSNG 

IF  (TYPE  .NE.  'C'  .OR.  LAYERS  .NE.  0 .OR.  WETHRS  . NE , 1)  GOTO  70 


READ  (INUNIT.  110)  N . Nl I . CL„ . H . CM . CH . W 
REREAD  120.  ( TEST ( I ) . 1=1 . 7 ) 


( TEST ( 1 ) 

.EO.  ' 

7 ) 

N-MISSNG 

(TEST (2) 

.EQ.  ' 

7 ) 

NH=MISSNG 

(TEST ( 3 ) 

.EQ.  ' 

7 ) 

CL.=MISSNG 

(TEST (4) 

.EQ.  ' 

7 ) 

H=M.TSSNG 

(TEST (5) 

.EQ.  ' 

7 ) 

CM®MISSNG 

< TEST ( 6 ) 

. EQ.  ' 

7 ) 

CH®MISSNG 

( TEST ( 7 ) 

.EQ.  7 7 

) 

W=MISSNG 

READ  ( INUNIT. 100»ERR=50)  TESTER 
REREAD  110.  UW 
GOTO  60 
50  UU=MISSNG 


60  READ  (INUNIT. 110)  PPP.  TT.  TDTD 


REREAD  120. 

( TEST  < I ) . ] 

[=1.3) 

IF  ( TEST ( 1 ) 

. EQ . 

7 7 ) 

PPP®  MI SSNG 

IF  < TEST  < 2 ) 

.EQ. 

7 7 ) 

TT®  MI SSNG 

IF  (TEST (3) 

.EQ. 

7 7 ) 

TDTD®  MI SSNG 

1B2 


CALL.  CONORT  < STATIN * JX  » JY  r JZ  > TIME  r I TIME ) 

WRITE  (OUTUNT >140)  JX » JY , JZ r IT IME r TYPEC t VISIBL r LAYERS 
WRITE:  ( OUTUNT  * 140 ) N»NH»CL»H»CM»CH»W 

WR I TE  < OUTUNT  r 140 ) WW , M I SSNG  rMISSNGfMl  SSNG  , M 1 SSNG . M I SSNG  • M 1 
GOTO  30 


70  TYPE  80 . DATE:  > T I ML  , ST  A I I N , TYPE » LAYERS » WETHRS  r 0 1 S I BL 
80  FORMAT  ('  ERROR  AT  ' »2I'7»2A5,3I4> 

90  CONTINUE 

100  FORMAT  (ID 
110  FORMAT  (01) 

120  FORMAT  (OA) 

130  FORMAT  (2A,3I> 

140  FORMAT  (8110) 

END 

SUBROUT I NE  CONOR!  ( ST AT  IN , JX » JY  r JZ . 1 1 ME r 1 I 1 ME ) 
INTEGER  JX  » JY  » JZ  i-  T I ME  t HOURS  r M I N 


HOURS® TIME/ 100 
MIN-TIME- (HOURS* 100) 
I r I ME- ( HOURS* AO ) EMI N 


IF  (STATIN  .NE 
JX-  5821 
JY®  2058 
JZ-  60 
GOTO  500 

10  IF  (STATIN  .NE 
JX-  7834 
JY®  4169 
JZ-  45 
GOTO  500 

20  IF  (STATIN  .NE 
JX-  6552 
JY-  4764 
JZ®  247 
GOTO  500 

30  IF  (STATIN  . NE 
JX®  5097 
JY®  46 
J Z-  1 1 

GOTO  500 

40  IF  (STATIN  .NE 

JX-  5527 
JY-  4400 
JZ®  315 
GOTO  500 

50  IF  (STATIN  .NE 
JX-  7531 
JY®  1801 


ABY ' ) GOTO  10 


AGS')  GOTO  20 


AHN')  GOTO  30 


AGO')  GOTO  40 


ATL')  GOTO  50 


AYS')  GOTO  60 


JZ*  46 
GOTO  500 

60  IF  (STATIN  .NE. 
JX*  7204 
JY  = 6 439 
JZ-  661 
GOTO  500 

70  II  (STATIN  .NE. 
JX-  3370 
JY-  4270 
JZ-  192 
GOTO  500 

00  IF  (STATIN  .NE. 
JX-  3431 
JY-  7106 
JZ-  104 
GOTO  500 

90  IF  (STATIN  .NE. 
JX-  016 
JY-  309 
JZ*  0 
GOTO  500 

100  IF  (STATIN  .NE. 
JX-  4708 
JY-  5919 
JZ*  210 
GOTO  500 

110  IF  (STATIN  .NE. 
JX-  7521 
JY - 5061 
JZ-  296 
GOTO  500 

120  IF  (STATIN  .NE. 
JX*  3360 
JY - 5479 
JZ-  196 
GOTO  500 

130  11  (STATIN  .NE. 
JX*  243 
JY*  2934 
JZ*  101 
GOTO  500 

140  11  (STATIN  .NE. 

JX*  6288 
JY*  3367 

jz*  no 

GOTO  500 

150  If  (STATIN  .NL. 
JX*  1499 
JY*  2924 
JZ*  94 


' AOL. ' ) GOTO  70 


KHfl  ) GOTO  80 


UNA  ) GOTO  90 


'HOE')  GOTO  100 


' CHA ' ) GOTO  110 


' GSP ' ) GOTO  120 


' MSV ' ) GOTO  130 


JAN  > GOTO  140 


MCN ' ) GOTO  150 


' hi.  1 > GOTO  160 


101 


CM  C4 


goto  500 

160  IF  (STATIN  .NE.  'MEM')  GOTO  170 
JX=  425 
JY-  5963 
JZ=  87 
GOTO  500 

170  IF  (STATIN  .NE.  ' MGM ' ) GOTO  180 
JX  = 3711 
JY-  2875 
JZ  = 62 

GOTO  500 

180  IF  (STATIN  . NE . 'MOB')  GOTO  190 
JX  = 1949 
JY^  1089 
JZ=  67 
GOTO  500 

190  IF  (STATIN  .NE.  ' MSY ' ) GOTO  200 
JX=  9 
JY-  351 
JZ==  9 
GOTO  500 

200  IF  (STATIN  .NE.  'F'NS')  GOTO  210 
JX-  2954 
JY=  842 
JZ-  36 
GOTO  500 

210  IF  (STATIN  .NE.  'SPA')  GOTO  220 
JX-  7764 
JY-  5891 
JZ=  251 
GOTO  500 

220  IF  (STATIN  . NE . ' TLH  ' ) GOTO  230 
JX=  5677 
JY - 779 

JZ=  21 
GOTO  500 

230  IF  (STATIN  .NE.  'TY5')  GOTO  240 
JX=  5873 
JY=  6815 
J Z=  299 
GOTO  500 

40  TYPE  250. STATIN 

50  FORMAT  ('  YOU  BLEW  IT  AT  STATION  '.AS) 

! 

500  CONTINUE 
RETURN 
END 


i h r» 
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PROGRAM  AIRWAY 

INTEGER  JX  f JY  f JZ  f TIME  f TYPE  f VISIBL  fLAYERS  fWETHRS  f I VISC  f ICLG  f ICLGV  r 
+ CO VAGE ( 1 0 ) f HE I GMT  ( 1 0 ) f TH I N ( 1 0 ) , WEATHR ( 7 ) f PRESSR f TEMP f 

+ DEWPT  f MISSNG f INUNIT  f OUTUNT 

LOGICAL  EOF 

DOUBLE  PRECISION  INFILE f OUTFIL 


COMMON  MISSNG  f I NUN IT  f OUTUNT 
MISSNG-  -32768 
EOF-  .FALSE. 

TYPE  10 

10  FORMAT  ('  NAME  THE  FILE. EXT  OF  THE  INPUT  FILE') 

ACCEPT  *.  INFILL 
TYPE  20 

20  FORMAT  ( ' TYPE  THE  LOGICAL  UNIT  NUMBER  OF  THE  INPUT  FILE') 
ACCEPT  *.  INUNIT 


TYPE  30 

30  FORMAT  ( ' TYPE  THE  OUTPUT  FILE. EXT') 
ACCEPT  *f  OUTFIL 


TYPE 40 

40  FORMAT  ('  TYPE  THE  LOGICAL.  UNIT  NUMBER  OF  THE  OUTPUT  FILE'  ) 
ACCEPT  *.  OUTUNT 

OPEN  (UNIT-  INUNIT f FILE-  INFILEf  ACCESS='SEQIN' > 

OPEN  ( UNIT-  OUTUNT.  FILE-  OUTFILf  ACCESS- ' SEQOUT ' ) 

50  CALL  ID  ( JX  f JY  f JZ  f TIME  f TYPE f VISIBL  f LAYERS f UETHRS f I VISC  f EOF  > 
IF  (EOF)  GOTO  60 

CALL  CLOUDS  ( LAYERS f ICLG f I CLOU f COO AGE f HEIGHT f THIN ) 

CALL  WHETHR  ( WEI  HRS . WEATHR f PRESSR f TEMP  t DEWPT  ) 

CALL  PR I NT  ( JX  f JY  f JZ  f T I ME . TYPE . V I S I BL f LAYERS f ICLG  f 
+ ICLGOf IOISCfWEATHRfCOVAGEfHEIGHTfTHIN) 

GOTO  50 
60  CONTINUE 
END 


SUBROUT I NE  ID  ( JX  f JY  f JZ  f T I ME  f TYPE  f V I S I BL  f L AYERSf WETHRS f I V I SC  f EOF ) 
INTEGER  JXf JYfJZf TIME f TYPE fUISIBL fLAYERS f WETHRS f I VISC f DATE 
REAL  VISIB 
LOGICAL  EOF 

DIMENSION  TEST (5) f LINE (20) 

COMMON  MISSNGf INUNI TfOUTUNT 

READ  ( INUNIT f40fEND=20)  ( TEST ( I ) f 1=1 f5 ) 

I F ( TEST ( A ) . EG . ' f ' ) GOTO  1 0 
IF  < TEST < 5 ) .EO.  'f')  GOTO  5 
TYPE  4 f < TEST ( I > f 1 = 1 f5 > 

4 FORMAT  ('  ERROR  'f5A2) 

DO  7 J=1 f 10 

READ  < INUNIT .9)  (LINE( I ) f 1=1 f20) 

TYPE  11  f ( LINE;  ( I ) f 1 = 1 f 20  ) 

7 CONTINUE 
STOP 

9 FORMAT  ( 20 A 1 ) 

11  FORMAT  ( IX f 20A1 ) 


1B6 


I ■ 


5 Kf.  Kt.AD  50  r DATE, TIME 
I IME- NEWT IM ( T IME ) 

READ  ( I NUN IT ,60)  < TEST ( 1 ) ,1-1,5) 

10  REREAD  60,  (TEST ( I ) , 1=1,5) 

REREAD  70,  STATIN , I TYPE , LAYERS , WETHRS , VISIB 

CALL  CONVRT  < STATIN , JX  , JY  , JZ  , ITYF'E  , TYPE  , OISIB  , VISIBL  , IVISC  ) 


IF 

(TEST (3) 

.EO. 

' ' ) 

LAYERS=MISSNG 

IF 

(TEST (4) 

.EQ. 

' ' ) 

WETHRS=MISSNG 

IF 

(TEST (5) 

.EO. 

' ' ) 

VISIBL  --MISSNG 

GOTO  30 
20  EOF- . TRUE . 

30  CONTINUE 
40  E OR MAT  (5A1) 

50  FORMAT  (2.1) 

60  FORMAT  ( 5 A ) 

7(  FORMAT  ( 2 A , 21 « F ) 

RE  TURN 
END 

INTEGER  FUNCTION  NEW  TIM  (TIME) 
I N r EGER  HOURS , M I N , T I ME 
HOUR S=T I ME/ 100 
MIN-TIME- ( HOURS# 100 ) 

NEWT I M= ( H0URS#60 ) +MI N 

RETURN 

END 


SUBROUTINE  CONVRT  ( STATIN , JX , JY  , JZ , ITYF'E , TYPE , OISIB , 
+ VISIBL, IVISC) 

INTEGER  JX , J Y , JZ , TYPE , T IME , I VI SC , HOURS , M I N , V I S I BL. 


REAL  VI SIB 

COMMON  MISSNG, INUNIT, OUTUNT 
DIMENSION  LINE (20) 

IF  (STATIN  « NE  . ' AE<Y  ' ) GOTO  10 
JX-  5821 
JY  = 2058 
JZ=  60 
GOTO  500 

10  IF  (STATIN  ,NE.  ' AGS ' ) GOTO  20 


JX=  7834 
JY=  4169 
JZ*  45 
GOTO  500 

20  IF  (STATIN  .NE.  ' AHN ' ) GOTO  30 
JX=  6552 
JY=  4764 
JZ-  247 
GOTO  500 

30  IF  (STATIN  . NE . 'AMG')  GOTO  40 
JX=  9314 
JY=  2207 
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Jl-  63 
GOTO  500 

40  IF  (STATIN  .NE.  ' ANB ' ) GOTO  50 
JX-  4214 
JY=  4302 
JZ  = 188 

GOTO  500 

50  IF  (STATIN  .NE.  'AND')  GOTO  60 
JX-  7082 
JY-  5396 
JZ  = 236 

GOTO  500 

60  IF  (STATIN  .NE.  ' AQQ ' ) GOTO  70 
JX=  5097 
JY-  46 
JZ  = 1 1 

GOTO  500 

70  IF  (STATIN  .NE.  ' ATI. ' ) GOTO  80 
JX=  5527 
JY  — 4400 
JZ=  315 
GOTO  500 

80  IF  (STATIN  .NE.  'AYS')  GOTO  90 
JX=  7531 
JY=  1801 
JZ=  46 
GOTO  500 

90  IF  (STATIN  . NE . ' AOL. ' ) GOTO  100 
JX=  7204 
JY=  6439 
JZ=  661 
GOTO  500 

100  IF  (STATIN  .NE.  ' BUM ' > GOTO  110 
JX-  3378 
JY=*  4278 
JZ=  .192 
GOTO  500 

110  IF  (STATIN  .NE.  ' BNA ' ) GOTO  120 
JX-  3431. 

JY=  7106 
JZ-  184 
GOTO  500 

120  IF  (STATIN  .NE.  'BTR'l  GOTO  130 
JX=  -837 
JY-  989 
JZ=*  23 
GOTO  500 

130  IF  (STATIN  .NE.  ' BOE ' ) GOTO  140 
JX-  816 
JY—  -389 
JZ-  0 


188 


GOTO  500 

140  IF  (STATIN  . NE . 
JX*  3609 
JY*  1194 

JZ*  56 
GOTO  500 

150  IF'  (STATIN  . NE . 
JX*  4708 
JY---  5919 
JZ*  210 
GOTO  500 

160  IF  (STATIN  .NE. 
JX*  2804 
JY*  7642 
JZ*  166 
GOTO  500 

170  IF  < ST At IN  .NE. 
JX*  5088 
JY*  3133 
JZ*  120 
GOTO  500 

180  IF  (STATIN  .NE. 
JX*  4875 
JY*  6938 
JZ*  570 
GOTO  500 

190  IF  (STATIN  .NE. 
JX*  4621 
JY*  1794 
JZ*  113 
GOTO  500 

200  IF  (STATIN  .NE. 
JX*  983 
JY*  7021 
JZ*  105 
GOTO  500 

210  IF  (STATIN  .NE. 
JX*  5446 
JY*  4546 
JZ*  257 
GOTO  500 

220  IF  (STATIN  . NE . 
JX*  556 
JY*  4256 
JZ*  40 


' CEW ' ) GOTO  150 


'CHA')  GOTO  160 


'CKO')  GOTO  170 


' CSG ' ) GOTO  180 


'CSO')  GOTO  190 


' DHN ' ) GOTO  200 


' DYR ' ) GOTO  210 


' F T Y ' ) GOTO  220 


' GLM ' ) GOTO  230 


GOTO  500 

230  IT  (STATIN  .NE.  'GNO'T  GOTO  240 
JX*  7731 
JY*  86 
JZ*  50 
GOTO  500 


lll«> 


GSP ' ) GOTO  250 


r 


240  IF  (STATIN  . NE . 
JX-  7521 
JY-  58  6 1 
,.J I - 296 

GOTO  500 

250  IF  (STATIN  . NE . 
JX  = 173 

J Y - 4250 
JZ-  41 
GOTO  500 

260  IF  (STATIN  .NE. 
JX-  3360 
J Y — 5479 
JZ=  196 
GOTO  500 

270  IF  (STATIN  . NE . 
JX-  243 
JY=  2934 


' GWO ' ) GOTO  260 


' HSV ' ) GOTO  270 


MAN')  GOTO  200 


JZ-  101 
GOTO  500 

200  IF  (STATIN  .NE. 
JX-  151 
J Y-  6053 
JZ-  805 
GOTO  500 

290  IF  (STATIN  . NE . 
JX-  -155 
JY-  1700 
JZ  = 143 

GOTO  500 

300  IF  (STATIN  .NE. 
JX-  6208 
JY=  3367 
JZ-  110 
GOTO  500 

310  IF  (STATIN  . NE . 
JX-  1499 
JY-  2924 
JZ-  94 
GOTO  500 

320  IF  (STATIN  .NE. 
JX-  425 
JY-  5963 
JZ=  07 
GOTO  500 

330  IF  (STATIN  .NE. 
JX-  3711 
JY - 2875 
JZ-  62 
GOTO  500 

340  IF  (STATIN  . NE . 


' JBR ' ) GOTO  290 


'MCE*')  GOTO  300 


'MCN' ) GOTO  310 


'MEI')  GOTO  320 


'MEH')  GOTO  330 


' MGM ' ) GOTO  340 


' MGFv ' ) GOTO  350 


I * 
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JX=  6200 
JY=  1569 
JZ=  88 
GOTO  500 

350  IF  (STATIN  .NE. 
JX=  1410 
JY  = 6549 
JZ=  129 
GOTO  500 

360  IF  (STATIN  .NE. 
JX-  1949 
JY=  1089 
JZ=  67 
GOTO  500 

370  IF  (STATIN  .NE. 
JX=  2582 
JY=  5592 
JZ=  171 
GOTO  500 

380  IF  (STATIN  .NE. 
JX=  9 
JY=  351 
JZ=  9 
GOTO  500 

390  IF  (STATIN  .NE. 
JX=  220 
JY=  400 
JZ=  3 
GOTO  500 

400  IF  (STATIN  .NE. 
JX=  2954 
JY=  842 
JI-  36 
GOTO  500 

410  IF  (STATIN  .NE. 
JX=  4833 
JY=  5162 
JZ=  196 
GOTO  500 

420  IF  (STATIN  .NE. 
JX=  8769 
JY=  -886 
JZ=  2 
GOTO  500 

430  IF  (STATIN  .NE. 
JX=  7764 
JY=  5891 
JZ=  251 
GOTO  500 

440  IF  (STATIN  .NE. 
JX=  2572 


' MKL ' ) GOTO  360 


'MOB')  GOTO  370 


' MSL ' ) GOTO  380 


' MSY ' ) GOTO  390 


' NEW ' ) GOTO  400 


' F’NS ' ) GOTO  410 


' RMG ' ) GOTO  420 


' SFB ' ) GOTO  430 


'SPA')  GOTO  440 


' TCL ' ) GOTO  450 


' 
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450 


4 AO 


4 A) 


480 


4 VO 


495 

497 


7 

9 

li 

too 


JY-  3910 
JZ=  57 
GOTO  500 
IF  (STATIN  • Nil  • 
JX-  9227 
JY-  1161 
JZ  = 1 1 

GOTO  500 
If  (STATIN  .NL'. 
JX  - 5677 
JY  = 77 9 

JZ  = 21 

GOTO  500 
IF  (STATIN  . NL. . 
JX-  1520 
JY-  5068 
JZ-  no 
GOTO  500 
IF  (STATIN  . NE . 
JX - 5873 
J Y-  6815 
JZ  = 299 

GOTO  500 
IF  (STATIN  • Nt . 


'TIX')  GOTO  460 


' 1 L H ' ) GOTO  4 70 


'TUP')  GOTO  400 


'TYS')  GOTO  490 


' OLD ' ) GOTO  495 


JX  = 6705 
JY-  1252 
JZ  = 66 

GOTO  500 
TYPE  497, STATIN 

FORMAT  ('  STATIN  ERROR  WITM',A5> 

DO  7 J-1,10 

READ  (I  NUN  IT,  9)  (LINE( I > , 1=1 » 20 > 
TYPE  11,  <LINE< I ) , 1=1 ,20) 
CONTINUE 
STOP 

FORMAT  ( 20A1 ) 


FORMAT  (IX , 20A1 ) 

IF  (I  TYPE  • E 0 « 'A')  TYPE=1 

IF  (I TYPE  . EO . 'S')  TYPE*  1 
V I S I BE  = I N T ( 0 1 S I D * 1 0 0 . ) 
IUISOMISSNG 
RETURN 


END 

SUBROUTINE  CLOUDS  (LAYERS , ICLG , ICLGV , COO AGE , HEIGHT • THIN) 
INTEGER  LAYERS , ICLG , ICLGV , COOAGE  ( 10 ) , HEIGHT < 10 ) , THIN<  10)  , L.AYEI 
DIMENSION  TEST (5) 

COMMON  MISSNG,  INUNIT, OLJTUNT 
ICLG=MISSNG 
ICLGO  = MISSNG 
DO  10  1=1,10 

THIN ( I ) =MISSNG 
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COVAGi:  < 1 ) =f1 1 SSNG 
HI  IGHT< 1 ) I SSNG 
10  CONTINUE 

U < L AYERS  .NE.  MISSNG)  GOTO  20 
Kt AH  < 1 NUN1 T f 50 ) < TEST < I ) , 1 = 1 »3) 

GO  10  40 

20  DO  40  LAYER- 1 r LAYERS 

READ  < INUNI T» 50)  < TEST < I ) r 1= 1 r 3 ) 

REREAD  60 f HEIGHT <1  AYER ) » COORGE  t HOWMES 
IF  ( TEST < 1 ) . EU.  ' ')  HEIGHT < LAYER ) -Ml SSNG 

IF  ( TEST < 2 ) .EU.  ' ')  COORGE-MISSNG 

IE  ( TEST ( 3 ) .EU.  ' ' .OR.  ICEG  .NE.  MISSNG)  GOTO  30 
ICLG=LAYER*10 

IF  (HOWMES  .EU.  'M'  .OR.  HOWMES  .EU.  'MO')  ICLG-ICLGf 1 

IF  (HOWMES  .EU.  'E'  .OR.  HOWMES  .EU.  'E0>  ICLG-ICLGF5 

IF  (HOWMES  .EU.  MO'  .OR.  HOWMES  . EO . 'EO')  ICLG0=1 


IF 

( CUORGE 

.EU. 

')  GOTO 

40 

IF 

(COORGE 

.EU. 

' CLR ' 

.OR. 

COORGE  .EU. 

' CLR 

) 

COOAGE  <L AYL  R ) -0 

ir 

< COORGE 

.EU. 

' SCT  ' 

.OR. 

COORGE  .EU. 

' SCI 

) 

CQOAGl  (1  AYE  R > 3 

IF 

i COORGE 

.EU. 

' DNN ' 

.OR. 

COORGE  .EU. 

' BI\N 

) 

COOAGE (LAYER)  - 6 

IF 

(COORGE 

.EU. 

'OOC' 

.OR. 

COORGE  .EU. 

' OOC  ' 

) 

COOAGE  U AYER ) -0 

IF 

( COORGE 

.EU. 

' -CLR ' 

.OR 

. COORGE  .EU. 

' SCT 

/ 

. UR . 

COORGE 

.EU. 

DRN  ' 

.OR 

. COORGE  .EU. 

' OOC 

) THIN(LAYER)-! 

40  CONTINUE 
50  FORMA)  ( 3A ) 

60  FORMAT  (I»2A) 

RE  TURN 
END 

SUBROUTINE  WHE I HR  ( WETHRS .WEATHR .PRESSR  » TEMP  . DEWI  I) 
INTEGER  WE  THRSf  UEATHR  ( 7 ) » F'RESSR  . TE  MR » DEWPT 
DIMENSION  TEST (3) 

COMMON  MISSNG. I NUN l T .OUTUNT 
DO  10  I-l»7 

WEA THR  < I ) -MISSNG 
10  CONTINUE 

IF  (WETHRS  .NE.  0 .AND.  WE  1 HRS  . NE . MISSNG) 

) READ  ( INUNIT. 20 > < WEATHR ( I ) . I = 1 . WETHRS ) 

READ  <INUN1T.30>  < TEST  Tl ) . I - 1 » 3 ) 

REREAD  40.  PRESSR  . TEMP  ► DEWF'T 


IF 

( TEST ( 1 ' 

.EU. 

' ' ) 

PRESSR=MISSNG 

IF 

( TEST (2) 

.EU. 

' ' ) 

TEMP-MISSNG 

IF 

< TE  ST  < 3 ) 

.EU. 

' ' ) 

DEWPT -MISSNG 

70  FORMAT  <71 ) 

30  FORMAT  < 3A ) 

40  FORMAT  (31) 

RETURN 
1 ND 

SUBROUTINE  PRINT  < JX . JY . JZ . T IME . TYPE , 0 I SIDE . LAYERS » I CEG » 

F ICLGO. IOISC. WEATHR. COOAGE. HEIGHT. THIN) 

INTEGER  JX. JY. JZ. TIME. TYPE. OISIPL .LAYERS. ICLGr ICLGO. OUTUNI . 
F IOISC. WEA THR < 7) . COOAGE < 10 ) . HE IGHT < 10 ) . THIN< 10 ) 


COMMON  M I SO  NO  * l NUN  l UOUTUNr 


UK  III  (Oil  II  IN  I.  10) 
UK  in  (OIJTUNI  , 10) 
WKIIl  v OU I UN  T » l 0 ) 
UK  III  (OU  1 1 IN  1.20) 
10  I ORMAT  (0110) 

20  FORMA  I (4  110) 

RE I URN 
l NH 


.IX  t JY  . J/  » 1 1 Ml  » r Yl’L  » V l S 1 BL  .LAY!  RO 

m o.  in  no,  i o i si: 

\WI  A 1 MR  ( l ) , 1-1 ,7) 

( LOOAUf  ( 1 ) » M 1 SONO , HEIGH T ( I ) , nil N ( I ) 
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